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ABSTRACT 

The purpose of this review is to describe in some detail the mathematical 
relationship between geometrodynamics and connection dynamics in the con- 
text of the classical theories of 2+1 and 3+1 gravity. I analyze the standard 
Einstein-Hilbert theory (in any spacetime dimension), the Palatini and Chern- 
Simons theories in 2+1 dimensions, and the Palatini and self-dual theories in 
3+1 dimensions. I also couple various matter fields to these theories and briefly 
describe a pure spin-connection formulation of 3+1 gravity. I derive the Euler- 
Lagrange equations of motion from an action principle and perform a Legendre 
transform to obtain a Hamiltonian formulation of each theory. Since constraints 
are present in all these theories, I construct constraint functions and analyze their 
Poisson bracket algebra. I demonstrate, whenever possible, equivalences between 
the theories. 
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1. Overview 

Einstein's theory of general relativity is by far the most attractive classical theory of 
gravity today. By describing the gravitational field in terms of the structure of spacetime, 
Einstein effectively equated the study of gravity with the study of geometry In general 
relativity, spacetime is a 4-dimensional manifold M with a Lorentz metric g a b whose curvature 
measures the strength of the gravitational field. Given a matter distribution described by 
a stress-energy tensor T a b, the curvature of the metric is determined by Einstein's equation 
G a b = 87r T ab . This equation completely describes the classical theory. 

As written, Einstein's equation is spacetime covariant. There is no preferred time vari- 
able, and, as such, no evolution. However, as we shall see in Section 2, general relativity 
admits a Hamiltonian formulation. The canonically conjugate variables consist of a positive- 
definite metric q a b and a density-weighted, symmetric, second-rank tensor field p ab — both 
defined on a 3-manifold E. These fields are not free, but satisfy certain constraint equations. 
Evolution is defined by a Hamiltonian, which (if we ignore boundary terms) is simply a sum 
of the constraints. 

Now it turns out that the time evolved data defines a solution, (M,g ao ), of the full 
field equations which is unique up to spacetime diffeomorphisms. In a solution, £ can be 
interpreted as a spacelike submanifold of M corresponding to an initial instant of time, while 
q ao and p ab are related to the induced metric and extrinsic curvature of £ in M.f\ Thus, the 
Hamiltonian formulation of general relativity can be thought of as describing the dynamics of 
3-geometries. Following Wheeler, I will use the phrase "geometrodynamics" when discussing 
general relativity in this form. 

On the other hand, all of the other basic interactions in physics — the strong, weak, and 
electromagnetic interactions — are described in terms of connection 1-forms. For example, 
the Hamiltonian formulation of Yang- Mills theory has a connection 1-form A a (which takes 
values in the Lie algebra of some gauge group G) as its basic configuration variable. The 
canonically conjugate momentum (or "electric field") E a is a density- weighted vector field 
which takes values in the dual to the Lie algebra of G. As in general relativity, these variables 
are not free, but satisfy constraint equations: The Gauss constraint D a E a = (where D a 
is the generalized derivative operator associated with A a ) tells us to restrict attention to 
divergence-free electric fields. Thus, just as we can think of the Hamiltonian formulation of 
general relativity as describing the dynamics of 3-geometries, we can think of the Hamiltonian 
formulation of Yang-Mills theory as describing the dynamics of connection 1-forms. I will 



2 More precisely, q a b is the induced metric on S, while p ab is related to the extrinsic curvature K a b via 
P ab = y/q(K ab - Kq ab ). 
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often use the phrase "connection dynamics" when discussing Yang-Mills theory in this form. 

Despite the apparent differences between geometrodynamics and connection dynamics, 
many researchers have tried to recast the theory of general relativity in terms of a connection 
1-form. Afterall, if the strong, weak, and electromagnetic interactions admit a connection 
dynamic description, why shouldn't gravity? Early attempts in this direction used Yang- 
Mills type actions, but these actions gave rise, however, to new theories of gravity A 
connection dynamic theory was gained, but Einstein's theory of general relativity was lost 
in the process. Later attempts (like the ones I will concentrate on in this review) left general 
relativity alone, but tried to reinterpret Einstein's equation in terms of the dynamics of a 
connection 1-form. The most familiar of these approaches is due to Palatini who rewrote 
the standard Einstein-Hilbert action (which is a functional of just the spacetime metric g a b) 
in such a way that the spacetime metric and an arbitrary Lorentz connection 1-form are 
independent basic variables. However, as we shall see in Section 6, the 3+1 Palatini theory 
does not succeed in recasting general relativity as a connection dynamical theory. The 3+1 
Palatini theory collapses back to the standard geometrical description of general relativity 
when one writes it in Hamiltonian form. 

More recently, Ashtekar [I], [| [| has proposed a reformulation of general relativity in 
which a real (densitized) triad Ef and a connection 1-form A l a (which takes values in the 
complexified Lie algebra of SO (3)) are the basic canonical variables. He obtained these 
new variables for the real theory by performing a canonical transformation on the standard 
phase space of real general relativity. For the complex theory, Jacobson and Smolin M 
and Samuel |5| independently found a covariant action that yields Ashtekar's new variables 
when one performs a 3+1 decomposition. This action is the Palatini action for complex 
general relativity viewed as a functional of a complex co-tetrad and a self-dual connection 
l-form.0 In one sense, it is somewhat surprising that these new variables could capture the 
full content of Einstein's equation since they involve only half (i.e., the self-dual part) of a 
Lorentz connection 1-form. On the other hand, the special role that self-dual fields play in 
the theory of general relativity was already evident in the work of Newman, Penrose, and 
Plebanski on self-dual solutions to Einstein's equation. In fact, much of this earlier work 
provided the motivation for Ashtekar's search for the new variables. 

Not only did the new variables give general relativity a connection dynamic descrip- 
tion; they also simplified the field equations of the theory — particularly the constraints. In 
terms of the standard geometrodynamical variables (q a b,P ab ), the constraint equations are 
non-polynomial. However, in terms of the new variables, the constraint equations become 

3 To recover the phase space variables for the real theory, one must impose reality conditions to select a 
real section of the complex phase space. 



3 



polynomial. This result has rekindled interest in the canonical quantization program for 3+1 
gravity. Due to the simplicity of the constraint equations in terms of these new variables, 
Jacobson, Rovelli, and Smolin || [7| and a number of other researchers have been able to 
solve the quantum constraints exactly. Although the quantization program has not yet been 
completed, the above results constitute promising first steps in that direction. 

The Palatini and self-dual theories described above were attempts to give general rel- 
ativity in 3+1 dimensions a connection dynamic description. A few years later, Witten 
H considered the 2+1 theory of gravity. He was able to show that this theory simplifies 
considerably when expressed in Palatini form. In fact, Witten demonstrated that the 2+1 
Palatini theory for vacuum 2+1 gravity was equivalent to Chern-Simons theory based on the 
inhomogeneous Lie group ISO(2, 1).[] He then used this fact to quantize the theory. This 
result startled both relativists and field theorists alike: relativists, since the Wheeler-DeWitt 
equation in geometrodynamics is as hard to solve in 2+1 dimensions as it is in 3+1 dimen- 
sions; field theorists, since a simple power counting argument had shown that perturbation 
theory for 2+1 gravity around a flat background metric is non-renormalizable — just as it is 
for the 3+1 theory. The success of canonical quantization and failure of perturbation theory 
in 2+1 dimensions came as a welcome surprise. Despite key differences between 2+1 and 
3+1 gravity (in particular, the lack of local degrees of freedom for 2+1 vacuum solutions), 
Witten's result has proven to be useful to non-perturbative approaches to 3+1 quantum 
gravity. In particular, since the overall structure of 2+1 and 3+1 gravity are the same (e.g., 
they are both diffeomorphism invariant theories, there is no background time, and the dy- 
namics is generated in both cases by 1st class constraints), researchers have been able to use 
2+1 gravity as a "toy model" for the 3+1 theory ||. 

Finally, the most recent developments relating geometrodynamics and connection dy- 
namics involve formulations of general relativity that are independent of any metric variable. 
This idea for 3+1 gravity dates back to Plebanski |I0||, and was recently developed fully 



by Capovilla, Dell, and Jacobson (CD J) ^ |E| [1J. Shortly thereafter, Peldan pf 
provided a similar formulation for 2+1 gravity. These pure spin-connection formulations 
of general relativity are defined by actions that do not involve the spacetime metric g a b in 
any way whatsoever — the action for the complex 3+1 theory depends only on a connection 
1-form (which takes values in the complexified Lie algebra of 50(3)) and a scalar density of 
weight —1. Moreover, the Hamiltonian formulation of this theory is the same as that of the 
self-dual theory, and by using their approach, CDJ have been able to write down the most 



Chern-Simons theory, like Yang-Mills theory, is a theory of a connection 1-form. However, unlike Yang- 
Mills theory, it is defined only in odd dimensions and does not require the introduction of a spacetime 
metric. 
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general solution to the 4 diffeomorphism constraint equations. Whether or not these results 
will lead to new insights for the quantization of the 3+1 theory remains to be seen. 

With this brief history of geometrodynamics and connection dynamics as background, 
the purpose of this review can be stated as follows: It is to describe in detail the theories 
mentioned above, and, in the process, clarify the mathematical relationship between ge- 
ometrodynamics and connection dynamics in the context of the classical theories of 2+1 and 
3+1 gravity. While preparing the text, I made a conscious effort to make the presentation as 
self-contained and internally consistent as possible. The calculations are somewhat technical 
and rather detailed, but I have included many footnotes, parenthetical remarks, and math- 
ematical digressions to fill various gaps. I felt that this style of presentation (as opposed to 
relegating the necessary mathematics to appendices at the end of the paper) was more in 
keeping with the natural interplay between mathematics and physics that occurs when one 
works on an actual research problem. Also, I felt that the added details would be of value 
to anyone interested in working in this area. 

In Section 2, 1 recall the standard Einstein-Hilbert theory and take some time to introduce 
the notation and mathematical techniques that I will use repeatedly throughout the text. 
Although this section is a review of fairly standard material, readers are encouraged to at least 
skim through the pages to acquaint themselves with my style of presentation. In Sections 3 
and 4, 1 restrict attention to 2+1 dimensions and describe the 2+1 Palatini and Chern-Simons 
theories and demonstrate the relationship between them. In Section 5, 1 couple a cosmo logical 
constant and a massless scalar field to the 2+1 Palatini theory. 2+1 Palatini theory coupled 
to a cosmological constant A is of interest since we shall see that the equivalence between 
the 2+1 Palatini and Chern-Simons theories continues to hold even if A ^ 0; 2+1 Palatini 
theory coupled to a massless scalar field is of interest since it is the dimensional reduction 
of 3+1 vacuum general relativity with a spacelike, hypersurface-orthogonal Killing vector 
field (see, e.g., Chapter 16 of [0). In fact, recent work in progress (by Ashtekar and 
Varadarajan) in the hamiltonian formulation of this reduced theory indicates that its non- 
perturbative quantization is likely to be successful. In Sections 6 and 7, I turn my attention 
to 3+1 dimensions and describe the 3+1 Palatini and self-dual theories. In Section 8, I 
couple a cosmological constant and a Yang-Mills field to 3+1 gravity. Section 9 describes 
a pure spin-connection formulation of 3+1 gravity, and Section 10 concludes with a brief 
summary and discussion of the results. All of the above theories are specified by an action. I 
obtain the Euler-Lagrange equations of motion by varying the action and perform a Legendre 
transform to put each theory in Hamiltonian form. I emphasize the similarities, differences, 
and equivalences of the various theories whenever possible. While this paper is primarily a 
review, some of the material is in fact new, or at least has not appeared in the literature in 
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the form given here. Much of Sections 3, 4, and 5 on the 2+1 theory fall in this category. 

I should also list a few of the topics that are not covered in this review. First, I have 
restricted attention to the more "standard" theories of 2+1 and 3+1 gravity. I have made 
no attempt to treat higher- derivative theories of gravity, supersymmetric theories, or any of 
their equivalents. Second, I have chosen to omit any discussion of quantum theory, although 
it is here, in quantum theory, that the change in emphasis from geometrodynamics to con- 
nection dynamics has had the most success. All of the theories described in this paper are 
treated at a purely classical level; issues related, for instance, to quantum cosmology and the 
non-perturbative canonical quantization program for 3+1 gravity are not dealt with. This re- 
view serves, instead, as a thorough pre-requisite for addressing the above issues. Moreover, 
many books and review articles already exist which discuss the quantum theory in great 
detail. Interested readers should see, in addition to the text books 0, [|, review articles 



17, 18, IS, PQ] and references mentioned therein. Third, in 2+1 dimensions, I have chosen 



to concentrate on the relationship between the 2+1 Palatini and Chern-Simons theories, and 
have all but ignored the equally interesting relationships between these formulations and the 
standard 2+1 dimensional Einstein-Hilbert theory. Fortunately, other researchers have al- 
ready addressed these issues, so interested readers can find details in |^T], |22|, Also, since 
Chern-Simons theory is not available in 3+1 dimensions, the equivalence of the 2+1 Palatini 
and Chern-Simons theories does not have a direct 3+1 dimensional analog. However, recent 



work by Carlip [^4], ^IJ and Anderson on the problem of time in 2+1 quantum gravity 
may shed some light on the corresponding issue facing the 3+1 theory. Finally, Section 9 on 
general relativity without-the metric deals exclusively with 3+1 gravity. Readers interested 
in a pure spin-connection formulation of 2+1 gravity should see |T5| . 

Penrose's abstract index notation will be used throughout. Spacetime and spatial tensor 
indices are denoted by latin letters from the beginning of the alphabet a,b,c, ■ ■ ■ , while 
internal indices are denoted by latin letters from the middle of the alphabet i,j, k, ■ ■ ■ or 
I, J,K, ■ ■ ■ . The signature of the spacetime metric g a b is taken to be ( — h +) or ( — h ++), 
depending on whether we are working in 2+1 or 3+1 dimensions. If V a denotes the unique, 
torsion- free spacetime derivative operator compatible with the spacetime metric g a b, then 
R a bc d kd '■= 2V[aVb]k c , R a b := R ac b c , and R := R a b9 ab define the Riemann tensor, Ricci tensor, 
and scalar curvature of V a . 

Finally, since I eventually want to obtain a Hamiltonian formulation for each theory, I 
will assume from the beginning that the spacetime manifold M is topologically £ x R. If the 
theory depends on a spacetime metric, I assume £ to be spacelike; if the theory does not 
depend on a spacetime metric, I assume S to be any (co-dimension 1) submanifold of M. In 
either case, I ignore all surface integrals and avoid any discussion of boundary conditions. 
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In this sense, the results I obtain are rigorous only for the case when £ is compact. Readers 
interested in a detailed discussion of the technically more difficult asymptotically flat case 
(in the context of the standard Einstein-Hilbert or self-dual theories) should see Chapters 
II.2 and III.2 of §. 

2. Einstein-Hilbert theory 

In this section, we will describe the standard Einstein-Hilbert theory. We obtain the vac- 
uum Einstein's equation starting from an action principle and perform a Legendre transform 
to put the theory in Hamiltonian form. We shall see that the phase space variables con- 
sist of a positive-definite metric q a b and a density-weighted, symmetric, second-rank tensor 
field p ab . These are the standard geometrodynamical variables of general relativity. We will 
also analyze the motions on phase space generated by the constraint functions and evaluate 
their Poisson bracket algebra. This section is basically a review of standard material. Our 
treatment will follow that given, for example, in Appendix E of |27J or Chapter II. 2 of 0. 

The standard Einstein-Hilbert theory is, of course, valid in n+1 dimensions. Everything 
we do in this section will be independent of the dimension of the spacetime manifold M. 
This is an important feature which will allow us to compare the standard Einstein-Hilbert 
theory with the Chern-Simons and self-dual theories. Unlike the standard Einstein-Hilbert 
theory, Chern-Simons theory is defined only in odd dimensions, while the self-dual theory is 
defined only in 3+1 dimensions. 

2.1 Euler-Lagrange equations of motion 

Let us begin with the well-known Einstein-Hilbert action 

S E H(g ab ) ■= I V=gR- (2.1) 
Jm 

Here g denotes the determinant of the covariant metric g a b, and R denotes the scalar cur- 
vature of the unique, torsion-free spacetime derivative operator V a compatible with g ab . I 
have taken the basic variable to be the contravariant spacetime metric g ab for convenience 
when performing variations of the action. The Einstein-Hilbert action is second- order since 
R contains second derivatives of g a b- 

To obtain the Euler-Lagrange equations of motion, we vary the action with respect to 
the field variable g ab . If we write the integrand as ^J^gR a i,g ab and use the fact that Sg = 
-9 9ab5g ab , we get 

SS EH = f V=g(Ra b - \Rga b )8g ab + / ^5R ab g ab . (2.2) 

JM 2 JM 
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The first integral is of the desired form, while the second integral requires us to evaluate the 
variation of the Ricci tensor R ab . Since one can show that^j 



SR ab g ab = V a v a (2.3) 

(where v a = ^ a {giJ>g bc ) — Vfe5g afc ), we see that modulo a surface integral, SSeh = if and 
only if 

G ab :=R ab - l -Rg ab = 0. (2.4) 

This is the desired result: The vacuum Einstein's equation can be obtained starting from 
an action principle. 

I should note that, strictly speaking, the variation of (2.1) with respect to g ab does not 
yield the vacuum Einstein's equation G ab = 0. The surface integral does not vanish since 
v a involves derivatives of the variation 5g ab . Even though Sg ab is required to vanish on the 
boundary, these derivatives need not vanish. This seems to pose a potential problem, but it 
can handled by simply adding to (2.1) a boundary term which will (upon variation) exactly 
cancel the surface integral. As shown in Appendix E of P7[, this boundary term involves the 



trace of the extrinsic curvature of the boundary of M. For the sake of simplicity, however, 
we will continue to use the unmodified Einstein-Hilbert action (2.1) and ignore all surface 
integrals as mentioned at the end of Section 1. 



2.2 Legendre transform 

To put the standard Einstein-Hilbert theory in Hamiltonian form, we will follow the 
usual procedure: We assume that M = E x R for some spacelike submanifold E and assume 
that there exists a time function t (with nowhere vanishing gradient (dt) a ) such that each 
t = const surface E t is diffeomorphic to E. To talk about evolution from one t = const surface 
to the next, we introduce a future-pointing timelike vector field t a satisfying t a (dt) a = 1. t a 
is the "time flow" vector field that defines the same point in space at different instants of 
time. We will treat t a and the foliation of M by the t = const surfaces as kinematical (i.e., 
non-dynamical) structure. Evolution will be given by the Lie derivative with respect to t a . 

Since we have a spacetime metric g ab as one of our field variables, we can also introduce a 
unit covariant normal n a and its associated future-pointing timelike vector field n a = g ab n b . 

5 To obtain this result, consider a 1-parameter family of spacetime metrics g a fc(A) and their associated 
spacetime derivative operators A V a . Define C ab c by A V a fc& =: V a fcf, + \C ab c k c and differentiate A V a <?f, c (A) = 
with respect to A. Evaluating this expression at A = gives C ab c = —^g cd C^aSg b d + ^ b 5g a d — ^d,Sg ab ), where 
g ab ■= g ab (0) and 5g ab := ^| A=0 ffah(A). Since R abc d W = R abc d + A 2W [a C b]c d + A 2 [C a ,C b ] c d , it follows 
that 5R ac := ix| A=0 #ab c b (A) = 2V [a C b]c b . Contracting with g ac (using Sg ab = -g ac g bd 5g cd ) yields the above 
result. 
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Note that since n a n a = — 1, g£ := 5 b + n a n b is a projection operator into the t = const 
surfaces. We will construct the configuration variables associated with the field variable g ab 
by contracting with n a and q b . We define the induced metric q ab , the lapse N, and shift N a 
via 

Qab ■■= qTQb 9mn (= 9ab + n a n b ) , 

N := -n a t b g ab , and (2.5) 
N a := ql t b . 

Note that in terms of N and iV a , we can write t a = Nn a + iV a . Furthermore, since N a n a = 
and q a bn a = 0, N a and q ab are (in 1-1 correspondence with) tensor fields defined intrinsically 
on S. 

The next step in constructing a Hamiltonian formulation of the Einstein-Hilbert theory 
is to decompose the Einstein-Hilbert action and write it in the form 

S EH {g ab )= JdtL EH (q,q). (2.6) 

Leh will be the Einstein-Hilbert Lagrangian provided it depends only on (q ab , N, N a ) and 
their first time derivatives. But, as written, (2.1) is not convenient for such a decomposition. 
The integrand yj '—gR contains second time derivatives of the configuration variable q ab . 
However, as we will now show, these terms can be removed from the integrand by subtracting 
off a total divergence. 

To see this, let us write the scalar curvature R as R = 2(G ab — R ab )n a n b . Then the 
differential geometric identities 

G ab n a n b = \{H- K ab K ab + K 2 ) and 

2 (2.7) 

R ab n a n b = -K ab K ab + K 2 + V b (n a V a n b - n b V a n a ) 

(where K ab := q™q b V m n n is the extrinsic curvature of the t = const surfaces and 1Z is the 
scalar curvature of the unique, torsion-free spatial derivative operator D a compatible with 
the induced metric q ab ) imply 

R = (K + K ab K ab - K 2 ) + (total divergence term). (2.8) 

Using the fact that yj— g = N^fq dt (where q denotes the determinant of q ab ), the Einstein- 
Hilbert action (2.1) becomes 

S EH (g ab ) = J dt y/qN(Jl + K ab K ab - K 2 ) + (surface integral). (2.9) 

If we ignore the surface integral, we get 

Leh = / s VqN(K + K ab K ab - K 2 ). (2.10) 
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This is the desired Einstein-Hilbert Lagrangian first proposed by Arnowitt, Deser, and Misner 
(ADM) [p8fl . The identity K a b = jkj(£t lab — 2D( a Nb)) allows us to express Leh in terms of 
only {q a bi N, N a ) and their first time derivatives. 

Given the Einstein-Hilbert Lagrangian, we are now ready to perform the Legendre trans- 
form. But before we do this, it is probably worthwhile to make a detour and first review the 
standard Dirac constraint analysis for a theory with constraints and recall some basic ideas 
of symplectic geometry. I propose to examine, in detail, a simple finite-dimensional system 
described by a Lagrangian 

L (q, q) ■= \q\ + 9392 - 94/(92, 93)- (2.11) 

Here (91, • • • , 94) G C are the configuration variables and (91, • • • , 94) are their associated time 
derivatives (or velocities). f(q2, 93) can be any (smooth) real-valued function of (52, 93). The 
techniques that arise when analyzing this simple system will apply not only to the standard 
Einstein-Hilbert theory but to many other constrained theories as well. Readers interested in 
a more detailed description of the general Dirac constraint analysis and symplectic geometry 
should see pj| and Appendix B of [Q, respectively. Readers already familiar with the 



standard Dirac constraint analysis may skip to the paragraph immediately following equation 
(2.25). 

To perform the Legendre transform for our simple system, we first define the momentum 
variables (pi, • • • ,04) via 

5L 

p 4 :=— (i=l, •••,4). (2.12) 

For the special form of the Lagrangian given above, they become 

Pi = 9i, P2 = 93, P3 = 0, and p A = 0. (2.13) 

Since only the first equation can be inverted to give q\ as a function of (q,p), there are 
constraints: Not all points in the phase space T = T*C = {(9i,Pi)| i = 1, •••,4} are 
accessible to the system. Only those (q, p) G r which satisfy 

<fa. ■= P2 ~ 93 = 0, (j) 2 :=P3 = 0, and 3 :=p 4 = O (2.14) 

are physically allowed. The 0j's are called primary constraints and the vanishing of these 
functions define a constraint surface in Tq. It is the presence of these constraints that 
complicates the standard Legendre transform. 

Following the Dirac constraint analysis, we now must now write down a Hamiltonian for 
the theory. But due to (2.14), the Hamiltonian will not be unique. The usual definition 
Ho(q,p) '■= Y^,=\PiAi — L(q, q) does not work, since there exist q^s which cannot be written 
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as functions of q and p. If, however, we restrict ourselves to the constraint surface defined 
by (2.14), we have 

H (q,p) = ^Pi 2 + qJ(q2,q3). (2.15) 

Since the right hand side of (2.15) makes sense on all of T , H (q,p) actually defines one 
possible choice of Hamiltonian. However, as we will show below, this Hamiltonian is definitely 
not the only one. 

For suppose Ai, A 2 , and A 3 are three arbitrary functions on r . Then 

H T (q,p) : = H (q,p) + A x 0i + A 2 2 + A 3 </> 3 

1 . ( 2 - 16 ) 
= gPi + 94/(92, 93) + Ai(p 2 - qs) + A 2 p 3 + A 3 £>4 

is another function (defined on all of T ) that agrees with H (q,p) on the constraint surface. 
Hx(q,p) is called the total Hamiltonian, and it differs from H (q,p) by terms that vanish on 
the constraint surface. This non-uniqueness of the total Hamiltonian exists for any theory 
that has constraints. 

Given Hx{q,p), the next step in the Dirac constraint analysis is to require that the 
primary constraints (2.14) be preserved under time evolution — i.e., that 

i:=feffr}o»0 (1=1,2,3). (2.17) 

Here ~ means equality on the constraint surface defined by (2.14) and { , }o denotes the 
Poisson bracket defined by the natural symplectic structur^ 

Q = dpi A dqi + dp 2 A dq 2 + dp 3 A dq 3 + dp^ A dq± (2-18) 

on r . Equation (2.17) is equivalent to the requirement that the evolution of the system 
take place on the constraint surface. 

Evaluating (2.17) for the primary constraints (2.14), we find that {4>3,Ht}o ~ implies 

04:=/(g 2 ,g 3 )~O. (2.19) 

6 A symplectic manifold (or phase space) consists of a pair (ro,Oo), where Tq is an even dimensional 
manifold and Qq is a closed and non-degenerate 2-form. (i.e., oIQq = and Qq(v,w) = for all w implies 
v = 0.) Qo is called the symplectic structure and it allows us to define Hamiltonian vector fields and 
Poisson brackets: Given any real-valued function / : To — * R, the Hamiltonian vector field Xf is defined by 
— ixj^o := df. Given any two real-valued functions /, g : Tq — > R, the Poisson bracket {f,g}o is defined 
by {/, g}o :— —Q,(Xt,Xg) = —Xf(g). As a special case, if Tq = T*Co is the cotangent bundle over some 
71-dimensional configuration space Cg , then f2 = dpi A dqi + • • • + dp n A dq n is the natural symplectic structure 
on r associated with the chart (q,p). It follows that {/, g} — J27=i(§^ §jt ~ Jp" af^")' w hich is the standard 
textbook expression for the Poisson bracket of / and g. 
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The other Poisson brackets yield conditions on Ai and A 2 . 04 is called a secondary constraint, 
and for consistency we must also require that 



04 := {0 4 ,#t}o ~0. (2.20) 

Here « now means equality on the constraint surface defined by (2.14) and (2.19). Since 
one can show that (2.20) follows from the earlier conditions on Ai and A 2 , (2.14) and (2.19) 
constitute all the constraints of the theory. 

The final step in the Dirac constraint analysis is to take all the constraints 1; • • • , 4 
and evaluate their Poisson brackets. If a constraint 0j satisfies {0j, 4>j}o ~ for all (j)j, then 
4>i is said to be 1st class. If, however, {0i, 0j}o ^ for some <pj, then 0j and <pj are said 
to form a 2nd class pair. (In terms of symplectic structures and Hamiltonian vector fields, 
a constraint 0£ is 1st class with respect to the symplectic structure Qq if and only if the 
Hamiltonian vector field defined by Qq is tangent to the constraint surface defined by 
the vanishing of all the constraints.) The goal is to solve all the 2nd class constraints (and 
possibly some 1st class constraints) and obtain a new phase space (r, Q) where the remaining 
constraints (pulled-back to T) are all 1st class with respect to the Poisson bracket defined by 
Q. Evaluating {0j, <pj}o for our simple system, we find that 03 is the only 1st class constraint 
with respect to f2 . By solving the second class pair 0i = and 2 = 0, we get 



and 



= dpi A dq\ + dq 3 A dq 2 + dp± A dq± (2-21) 

!>i=0, 02=0 



En 



, nA =\pi 2 + q±fXq2,qz) + \m- (2.22) 

01=0, 02=0 Z 

The remaining constraints 03 and 04 are now both 1st class with respect to this new sym- 
plectic structure. 

Although we have successively eliminated all the 2nd class constraints, we can go one 
step further. We can solve the 1st class constraint 03 := p± = by gauge fixing the configu- 
ration variable q±. Even though this step is not required by the Dirac constraint analysis, it 
simplifies the final phase space structure somewhat. Solving 3 = and pulling-back (2.21) 
and (2.22) to this new constraint surface T (coordinatized by (qi, g 2 , <?3,Pi)), we obtain 

n := dp x A dqi + dq 3 A dq 2 (2.23) 

and 

H( qi ,q 2 ,q 3 , Pl ) := - Pl 2 + q 4 f(q 2 , q 3 ). (2.24) 

Here q± is no longer thought of as a dynamical variable — it is a Lagrange multiplier of the 
theory associated with the 1st class constraint f(q 2 ,q 3 ) = 0. 
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To summarize: Given a Lagrangian of the form 

L (q, q) ■= \q\ + 9392 - 94/(92, 93), (2.25) 

the Dirac constraint analysis says that the momentum p\ is unconstrained, while P2 = 93 
and P3 = Pa = 0. Demanding that the constraints be preserved under evolution, we obtain a 
secondary constraint f(q2, 93) = 0. The constraints p<i — q% = and P3 = form a 2nd class 
pair and are easily solved; the remaining constraints p^ = and f(q2, 93) = now form a 1st 
class set. By gauge fixing g 4 we can solve p^ = 0, and thus obtain a new phase space (T, Q) 
coordinatized by (91, 92, 93, £>i) with symplectic structure (2.23) and Hamiltonian (2.24). We 
are left with a single 1st class constraint, f(q2, 93) = 0. 

Let us now return to our analysis of the standard Einstein-Hilbert theory. Given Leh, 
we find that the momentum p ab canonically conjugate to q ab is given by 

pab ;= 5Leh_ = {Rab _ Rqab) ^ (2 2g) 

oCt q a b 

while the momenta canonically conjugate to N and iV a are zero. Since equation (2.26) can 
be inverted to give 

C r q ab = 2Nq- 1 ' 2 {p ab - ^pq ab ) + 2D (a N b) , (2.27) 

it does not define a constraint. However, iV and iV a play the role of Lagrange multipliers. 

Thus, by following the Dirac constraint analysis we find that the phase space (Teh, &eh) 
of the standard Einstein-Hilbert theory is coordinatized by the pair (q a b,p ab ) arid has sym- 
plectic structure^ 

= / dt^ fe A dq ab . (2.28) 



The Hamiltonian is given by 

H EH (q,p) = J E N {~ 1 1/2 K + q- 1,2 (p ab p ab - \p 2 )) - 2N a q ab D c f c . (2.29) 

As we shall see in the next subsection, this is just a sum of 1st class constraint functions 
associated with 

C(q,p) := -q 1/2 K + q~ 1/2 (p ab p ab - \f) « and (2.30a) 
C a (q,p) := -2q ab D c p hc w 0. (2.306) 



I use d and A to denote the infinite-dimensional exterior derivative and infinite-dimensional wedge 
product of forms on Teh- They are to be distinguished from d and A which are the finite-dimensional 
exterior derivative and finite-dimensional wedge product of forms on E. Note that in terms of the Poisson 
bracket { , } defined by VL E h, we have {q ab (x),p cd (y)} = 6?S^S(x,y). 
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Note that constraint equation (2.30a) is non-polynomial in the canonically conjugate vari- 
ables due to the dependence of 1Z on the inverse of q ab . This is a major stumbling block 
for the canonical quantization program in terms of (q a b,P ab )- To date, there exist no ex- 
act solutions to the quantum version of this constraint in full (i.e., non-truncated) general 
relativity. 

2.3 Constraint algebra 

To evaluate the Poisson brackets of the constraints and to determine the motions they 
generate on phase space, we must first construct constraint functions (i.e., mappings Teh — * 
R) associated with the constraint equations (2.30a) and (2.306). To do this, we smear C(q,p) 
and C a (q,p) with test fields N and N a on E — i.e., we define 

C(N) : = jf N( - q^K + q- 1/2 (F b Pa b - \f)) and (2.31a) 

C(N) := J -2N a q ab D c f? c . (2.316) 

They are called the scalar and vector constraint functions of the standard Einstein-Hilbert 
theory. 

The next step is to evaluate the functional derivatives of C(N) and C(N). For recall 
that if /, g : T E h —> R are any two real- valued functions on phase space, the Hamiltonian 
vector field Xf (defined by the symplectic structure (2.28)) is given by 

X = f 1LJ^_1LA_ (2 32) 

f 7s 5p ab 5q ab 5q ab Sp ab 

and the Poisson bracket {/, g} (defined by {/, g} := —Xf(g)) is given by 

U'V-hSqaSp* Sp^5q ab {2 - 66) 
Note that under the 1-parameter family of diffeomorphisms on Teh associated with Xf, 

q a b ^ qab + e—^ + 0(e 2 ) and (2.34a) 

p ab ^p ab -e^- + 0(e 2 ). (2.346) 
oq ab 

— * 

We will use (2.33) to determine the various Poisson brackets between C(N) and C(N); we 
will use (2.34a) and (2.346) to determine the motions that they generate on phase space. 

Let us begin with the vector constraint C(N). Integrating (2.316) by parts and noting 
that 2D(aN b ) = C^qab, we get 

C(N) = J^q ab )if b (= - jf q ab {C^ b )) ■ (2.35) 
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By inspection, 



Thus, we see that 



^-CJ? and (2,6) 



q a b >-> g a 6 + e£^g a b + 0(e 2 ) and (2.37a) 

~a6 ^ ~a6 + + Q(e 2 ) (2.376) 

is the motion on Teh generated by C(N). Note that (2.37a) and (2.376) are the maps on 
the tensor fields q ab and p ab induced by the 1-parameter family of diffeomorphisms on X 
associated with the vector field N a . In other words, the Hamiltonian vector field X c ,gs on 
Teh is the lift of the vector field N a on E. 

Let us now consider the scalar constraint C(N). Due to the non-polynomial dependence 
of 1Z on q ab , the functional derivative SC(N)/5q a b is much harder to evaluate. After a fairly 
long calculation, one finds that|] 

= -]-NC(q,p)q ab + 2iV<T 1 /2(pV c - h^) 
d qab 2 2 (2.38) 

+ Nq 1/2 (K ab - Kq ab ) - q 1/2 (D a D b N - q ab D c D c N). 
A much simpler calculation gives 

5 -^ = 2Nq-^{p ah - l -pq ab ). (2.39) 

Recall that for the vector constraint function C(N), the motion on along corre- 
sponded to the Lie derivative of g ab and p ab with respect to iV a . Thus, one might expect the 
motion on along Xq(n) to correspond to the Lie derivative with respect to t a := Nn a . 
We will now show that if we restrict ourselves to the constraint surface Teh C Teh (defined 
by (2.30a) and (2.306)), then this is actually the case. 

Comparing (2.39) with equation (2.27) (setting iV a = 0), we see that 5C(N)/5p ab = 

CtQab, SO 

q ab >-> q ab + eC r q ab + 0(e 2 ) (2.40) 

as conjectured. Similarly, writing p ab = ^/q(K ab — Kq ab ) and using the differential geometric 
identity 

C N aK ah = -N1l ab + 2NK a c K bc - NKK ab + D a D b N (2.41) 

8 To obtain this result we used the facts that Sq = q q ab 5q a b and 51Z a bq ab = D a v a for v a = —D a {q bc 8qb c ) + 
D b (q ac 5qbc)- These are just the spatial analogs of the results used in subsection 2.1 when we varied the 
Einstein-Hilbert action with respect to g ab . 
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(which holds in this form when R ab = 0), we see that 

5C(N) fl 



Sq, 



ab 



-(-NC(q,p)q ab + C^ b y (2.42) 



Thus, 

pab ^pab + e (^NC(q,p)q ab + C r p ab ^j + 0(e 2 ). (2.43) 
If we now restrict ourselves to T E h (so that C(q,p) = 0), we get 



This is the desired result: When restricted to Teh C Teh, the Hamiltonian vector field 
Xc(N) on Teh is the lift of the vector field t a := Nn a on E. 

We are now ready to evaluate the Poisson brackets between the constraint functions. But 
first note that if f(M) : T E h R is any real- valued function on phase space of the form 

f(M) := jf M a - b c ... d f a ... b c - d (q,P) (2.45) 

(were M a "' b c ... d is any tensor field on E which is independent of q ab and p ab ), then 

= jl .U" Cj a ... b c - d (q,p). 
Integrating the last line of (2.46) by parts and throwing away the surface integral, we get 

{C(N),f(M)} = f(CffM). (2.47) 

Thus, the Poisson bracket of C(N) with any other constraint function is easy to evaluate. 
We have 

{C{N),C(M)} = C{[N,M]) and (2.48a) 
{C(N),C(M)} = C(£$M), (2.486) 

— * — * 

where [N,M] := C^M a is the commutator of the vector fields iV a and M a on E. Note 
that (2.48a) tells us that the subset of vector constraint functions is closed under Poisson 
brackets. In fact, iV a i— > C(N) is a representation of the Lie algebra of vector fields on E. 
The commutator of vector fields on E is mapped to the Poisson bracket of the corresponding 
vector constraint functions. 



5q 



ab 



(2.46) 
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We are left with only the Poisson bracket {C(N), C(M)} of two scalar constraint func- 
tions to evaluate. Using (2.38) and (2.39) (and eliminating all terms symmetric in M and 
N), we get 



where K a := (Nd a M - Md a N) = q ab (Nd b M - Md b N). Thus, the Poisson bracket of two 
scalar constraints is a vector constraint. Although this implies that the subset of scalar 
constraint functions is not closed under Poisson bracket, the totality of constraint func- 
tions (scalar and vector) is — i.e., the constraint functions form a 1st class set as claimed in 
subsection 2.2. Note, however, that since the vector field K a depends on the phase space 
variable q ab (through its inverse), the Poisson bracket (2.49) involves structure functions. 
The constraint functions do not form a Lie algebra. 

3. 2+1 Palatini theory 

In this section, we will describe the 2+1 Palatini theory which, as we shall see at the end 
of subsection 3.2, is defined for any Lie group G. We will discuss the relationship between 
the Palatini and Einstein-Hilbert actions, and show how the 2+1 Palatini theory based 
on 5*0(2, 1) reproduces the standard results of 2+1 gravity. After performing a Legendre 
transform to put this theory in Hamiltonian form, we shall see that the phase space variables 
consist of a connection 1-form (which takes values in the Lie algebra of G) and its 
canonically conjugate momentum (or "electric field") Ef. Thus, for G = 50(2, 1), the 2+1 
Palatini theory gives us a connection dynamic description of 2+1 gravity. The constraint 
equations are polynomial in the basic variables and the constraint functions form a Lie algebra 
with respect to Poisson bracket. 

Once we write the 2+1 Palatini action in its generalized form, we will let G be an arbitrary 
Lie group. To reproduce the results of 2+1 gravity, we simply take G to be 50(2, 1). Note 
that much of the material in subsections 3.2 and 3.3 can also be found in ||30|| . 

3.1 Euler-Lagrange equations of motion 

Recall the standard Einstein-Hilbert action of Section 2, 



{C(N), C(M)} = j -2M(D a D b N - q ab D c D c N)(p ab - -pq ab ) - (N M) 
7s 2 




(2.49) 



C(K) 




(3.1) 
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To define the 2+1 Palatini action, it is convenient to first rewrite the integrand y/—gR in 
triad notation. But in order to do this, we will have to make a short mathematical digression. 
Readers interested in a more detailed discussion of what follows should see f3Tfl . Readers 
already familiar the method of orthonormal bases may skip to the paragraph immediately 
following equation (3.9). 

Consider an n-dimensional manifold M, and let V be a fixed n-dimensional vector space 
with Minkowski metric 7]u having signature ( — !-■••+). A soldering form at p e M is an 
isomorphism e^(p) : T p M — > V. (Here T p M denotes the tangent space to M at p.) Although 
an n-manifold does not in general admit a globally defined soldering form e„, we can use 
to define tensor fields locally on M. For instance, 

g a b ■■= elfiivu ( 3 - 2a ) 

is a (locally defined) spacetime metric having the same signature as rjjj. The inverse of 
will be denoted by ef; it satisfies 

gabe^j = Vjj . (3.2b) 

Spacetime tensor fields with additional internal indices I, J,K, - ■ ■ will be called generalized 
tensor fields on M. Spacetime indices are raised and lowered with the spacetime metric g a b] 
internal indices are raised and lowered with the Minkowski metric rjjj. 

If one introduces a standard basis {h\ \I_ — 1, ■ • • , n} in V, then the vector fields ef := e a jb\ 
form an orthonormal basis of g ab . These n-vector fields will be called a triad when n = 3 
and a tetrad when n = 4. The dual co- vector fields, e^ := g a bi]—e b j, will be called a co-triad 
and a co-tetrad when n = 3 and 4, respectively. I should note, however, that from now on 
I will ignore the distinction between a soldering form and the co-vector fields e^. I will 
call a 3-dimensional soldering form \ l a a co-triad and a 4-dimensional soldering form ^ a 
co-tetrad in what follows. 

To do calculus with these generalized tensor fields, it is necessary to extend the definition 
of spacetime derivative operators so that they also "act" on internal indices. We require 
(in addition to the usual properties that a spacetime derivative operator satisfies) that a 
generalized derivative operator obey the linearity, Leibnitz, and commutativity with contrac- 
tion rules with respect to the internal indices. Furthermore, we require that all generalized 
derivative operators be compatible with rjij. Given these properties, it is straightforward to 
show that the set of all generalized derivative operators has the structure of an affine space. 
In other words, if d a is some fiducial generalized derivative operator (which we treat as an 
origin in the space of generalized derivative operators) , then any other generalized derivative 
operator V a is completely characterized by a pair of generalized tensor fields A ab c and A aI J 



18 



defined by 

V a k bI =: d a k bI + A ab c k cI + A aI J k bJ . (3.3) 

We will call A ab and A aI J the spacetime connection 1-form and internal connection 1-form 
of V a . It is easy to show that 

A aIJ = A a[IJ] and A ab c = A {ab) c . (3.4) 

These conditions follow from the requirements that all generalized derivative operators be 
compatible with rju and that they be torsion-free. Later in this section, we will consider 
what happens if we allow derivative operators to have non-zero torsion — i.e., if A[ ab f ^ 0. 
Finally, note that A ab c need not equal A a j J e I b e c J , in general. 

As usual, given a generalized derivative operator V a , we can construct curvature tensors 
by commuting derivatives. The internal curvature tensor F abI J and the spacetime curvature 
tensor F abc d are defined by 

2V [a V b] kj =: F abI J kj and (3.5a) 
2V [a V b] k c =:F abc d k d . (3.56) 

If our fiducial generalized derivative operator is chosen to be flat on both spacetime and 
internal indices, then 

F abI J = 2d [a A b]I J +{A a ,A b } I J and (3.6a) 
F abc d = 2d [a A b]c d +[A a ,A b } c d . (3.66) 

Here [A a , A b ]/ := {A aI K A bK J - A bI K A aK J ) and [A a , A b ] d := {A ac e A be d - A bc e A ae d ) are the 
commutators of linear operators. 

Just as a compatibility with a spacetime metric g ab defines a unique, torsion-free spacetime 
derivative operator V a , compatibility with an orthonormal basis (and thus with g ab ) defines 
a unique torsion- free generalized derivative operator, which we also denote by V a . The 
Christoffel symbols T aI J and T ab c are defined by 

V a k bI =: d a k bI + T ab c k cI + r aI J k bJ , (3.7) 

and satisfy 

T a i J = -e bJ (d a e bI + T ab c e cI ) and (3.8a) 

T ab c = ~g cd (d a g bd + d b g ad - d d g ab ). (3.86) 
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It also follows that internal and spacetime curvature tensors R a bi J and R a bc d of V a are related 
by 

R a bi J = Rab/ejej. (3.9) 

We will need the above result in this and later sections to show that the Palatini and self-dual 
actions reproduce Einstein's equation. 

Now let us return to our discussion of the 2+1 Palatini theory. Recall that we wanted to 
write the integrand ^—gR in triad notation. Using 

R a bi J = R a bc d 3 e/ 3 e J d (3.10) 

(which is equation (3.9) written in terms of a triad 3 ej) and 

e abc = 3 el 3 ei 3 e« e IJK (3.11) 

(which relates the volume element e a bc of g a b = V a 3 e b Vu to the volume element euk of rju), 
we find that 

V^gR = V^g S b [d 5 c e] R bc de 

_ 1 Z^abc p de 

1 (3-12) 

^.abc 3 a I 3 a J 3 a K p de 

2^ a d e eiJKKbc 

— ^n ahc c V R JK 

— 2^ e UK e a Kbc 

Thus, viewed as a functional of a co-triad 3 e^, the standard Einstein-Hilbert action is given 
by 

SehCc) = \j jf c e IJK % R bc JK . (3.13) 

To obtain the 2+1 Palatini action, we simply replace R ab i J in (3.13) with the internal 
curvature tensor 3 F abI J of an arbitrary generalized derivative operator ^D a defined by 

^Dakj-.^dakj + ^/kj. (3.14) 

We define the 2+1 Palatini action based on SO(2, 1) to be 

S P { 3 e 3 A):=\f fj ab % JK ¥ a 3 F bc JK , (3.15) 

4 JM 

where 3 F ab i J = 2d[ a 3 A b ]i J + [ 3 A a , 3 A b ]i J . Note that I have included an additional factor of 
1/2 in definition (3.15). This overall factor will not affect the Euler-Lagrange equations of 
motion in any way, but it will change the canonically conjugate variables. I have chosen 
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to use this action so that the expressions for our canonically conjugate variables agree with 
those used in the literature (see, e.g., |30| ). 

As defined above, Sp( 3 e, 3 A) is a functional of both a co-triad V a and a connection 1-form 
3 A a j J which takes values in the defining representation of the Lie algebra of SO(2, 1). Note 
also that ^V> a as defined by (3.14) knows how to act only on internal indices. We do not require 
that ^Da know how to act on spacetime indices. However, when performing calculations, we 
will find that it is often convenient to consider a torsion-free extension of to spacetime 
tensor fields. It turns out that all calculations and all results will be independent of our 
choice of torsion- free extension. In fact, we will see that these results hold for extensions of 
that have non-zero torsion as well. 

Since the 2+1 Palatini action is a functional of both a co-triad and a connection 1-form, 
we will obtain two Euler-Lagrange equations of motion. When we vary V a1 we get 

V abc euK 3 F bc JK = 0. (3.16) 

When we vary 3 A a IJ , we get 

"Vbiv^euK 3 ef ) = 0. (3.17) 

To arrive at (3.17), we considered a torsion-free extension of ^Da to spacetime tensor fields 
(so that 5 3 F bc JK = 2 ^D^A^ JK ) and then integrated by parts. The surface integral vanished 
since 5 3 A C JK = on the boundary, while the remaining term gave (3.17). Note that since 
the left hand side of (3.17) is the divergence of a skew spacetime tensor density of weight +1 
on M, it is independent of the choice of torsion- free extension of Since fj abc ejjK 3 ef = 
2( 3 e) 3 e^ 3 e b j (where ( 3 e) := y/—g), we can rewrite (3.17) as 

*D h (( 3 e) 3 e' a 3 e?) = 0. (3.18) 

We shall see in Section 6 that the form of equation (3.18) holds for the 3+1 Palatini theory 
as well. 

To determine the content of equation (3.18), let us express in terms of the unique, 
torsion- free generalized derivative operator V a compatible with 3 e 7 , and ^a/" 7 defined by 

V a k I + 3 C aI J kj. (3.19) 

Since (3.18) is the divergence of a skew spacetime tensor density of weight +1 on M, and 
since V a is compatible with 3 e 7 , we get 

3 C b i K 3 * 3 ej + 3 C hJ K 3 ef 3 4 = 0. (3.20) 

This is equivalent to the statement that the (internal) commutator of ^Cmj and 3 ej 3 e b j 
vanishes. We will now show that (3.20) implies that 3C a / = 0. 
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To see this, define a spacetime tensor field 3 S a b c via 

Xbc := 3 C aI j 3 e J b \{. (3.21) 

(Note, incidently, that 3 S a b c is not the spacetime connection of relative to V .) Then the 
condition ^C a /j = ^C a [ij\ is equivalent to 3 S a bc = 3 S a [bc\- Now contract equation (3.20) with 
3 e I a 3 e J c . This yields 3 Sb c b = 0, so 3 S a b c is trace-free on its first and last indices. Using this 
result, (3.20) reduces to 

3 C bI K 3 e a K 3 e b j - 3 C bJ K 3 e b j 3 e a K = 0. (3.22) 
If we now contract (3.22) with 3 e^ 3 ej, we get 

bed = ^(cd) ■ {6.26) 

Thus, 3 S a bc is symmetric in its first two indices. Since 3 S a b c = 3 S a [b c ] and 3 S a b c = 3 S( a b) c , we 
can successively interchange the first two and last two indices (with the appropriate sign 
changes) to show 3 S a b c = 0. Futhermore, since e\ are invertible, we get 3 C a i J = 0. This is 
the desired result .[] 

Since 3 C a i J = 0, we can conclude that the generalized derivative operator must agree 
with V a when acting on internal indices. Thus, although the Palatini action started as a 
functional of a co-triad and an arbitrary generalized derivative operator ^Dq, we find that 
one equation of motion implies that = V a . In terms of connection 1-forms, ^q/* 7 = 
implies that 3 A a j J = T a j J , where T a j J is the internal Christoffel symbol of V a . Using this 
result, the remaining Euler-Lagrange equation of motion (3.16) becomes 

V abc eiJKR bc JK = 0. (3.24) 

When (3.24) is contracted with 3 e dI , we get G ad = 0. Thus, the Palatini action based on 
5*0(2, 1) reproduces the standard 2+1 vacuum Einstein's equation. 

It is interesting to note that to show that the Palatini action reduces to the standard 
Einstein-Hilbert action in 2+1 dimensions, we need only vary the connection 1-form 3 A a i J . 
Since we found that (3.17) could be solved uniquely for 3 A a i J in terms of the remaining basic 
variables 3 e^, we can pull-back Sp( 3 e, 3 A) to the solution space 3 A aI J = T aI J and obtain a new 
action S_p( 3 e), which depends only on a co-triad. This pulled-back action is just 1/2 times 
the standard Einstein-Hilbert action 5W( 3 e) given by (3.13). But what about the boundary 
term that one should strictly include in the standard Einstein-Hilbert action? It looks as if 
S_p{ 3 e) is missing this needed term. 

9 This method of proving 3 C a i J — — which generalizes to the 3+1 Palatini and self-dual actions — was 
shown to me by J. Samuel and A. Ashtekar. 
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The answer to this question is the following: Whereas the standard Einstein-Hilbert 
action is a second-order action, the 2+1 Palatini action is first-order. As mentioned at the 
beginning of Section 2, varying the standard Einstein-Hilbert action (3.1) with respect to g ab 
produces a surface integral involving derivatives of the variation Sg ab . Since we are allowed 
only to keep g ab fixed on the boundary, this surface integral is non-vanishing and must be 
compensated for by adding a boundary term to (3.1). This is also the case if we vary Seh( 3 z) 
given by (3.13) with respect to V a . On the other hand, when we vary the Palatini action 
(3.15) with respect to 3 A a j J , we hold 3 A a j J fixed on the boundary and 3 e^ fixed throughout. 
Then by solving (3.17) uniquely for 3 A a j J , we can pull-back Sp( 3 e, 3 A) to the solution space 
3 A aI J = T aI J . But now when we vary S_ P ( 3 e) with respect to 3 e J a which lie entirely in the 
solution space, fixing 3 e„ on the boundary also fixes certain derivatives of 3 e„ on the boundary. 
This is a reflection of the fact that the reduction procedure comes with a prescription on how 
to do variations. It is precisely the vanishing of these derivatives of <5 3 e^ which eliminates 
the need of a boundary term for S_ P ( 3 e). 

It is also interesting to note that we could obtain the same result ( 3 A aI J = T aI J ) by 
considering an extension of ^Da to spacetime tensor fields with non-zero torsion ^T^ . (Recall 
that if 3 A a b c denotes the spacetime connection 1-form of the extension of ^Da, then the torsion 
tensor ^T^ is defined by 2 3 P [a *D h] f =: ^DJ and satisfies ^T^ = 2 3 A [ab f.) By varying 
the 2+1 Palatini action (3.15) with respect to 3 A aI J , we would find 

2 ^VeJ - % b c 3 e[ = 0. (3.25) 

This is the field equation for 3 A aI J which holds for any extension of ^D a to spacetime tensor 
fields. If we restrict ourselves to torsion- free extensions, we get back equation (3.17). Then 
by following the argument given there, we would find 3 A aI J = T aI J as before. 

However, there exists an alternative approach to solving equation (3.25) which is often 
used by particle physicists. Namely, instead of considering a torsion-free extension of to 
spacetime tensor fields, one considers an extension of to spacetime tensor fields which is 
compatible with the co-triad 3 e^. This can always be done, but the price of such an extension 
is in general a non-zero torsion tensor 3 T a b c . But since we now have 3 e{ = 0, equation 
(3.25) implies 

Xi, 3 e{ = 0. (3.26) 

Invertibility of 3 e\ then implies that 3 T a 5 c = 0. Since there exists only one torsion-free 
derivative operator compatible with 3 e^, we can conclude that = V a (or equivalently, 
3 A aI J = T aI J ). This is the desired result. 

Finally, to conclude this section, let us write the 2+1 Palatini action (3.15) in a form 
which is valid for any Lie group G. Recall that the connection 1-form 3 A aI J — being a linear 
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operator on the internal 3- dimensional vector space (equipped with the Minkowski metric 
rju) and satisfying 3 A a u = 3 A a [ij] — takes values in the defining representation of the Lie 
algebra of 5*0(2,1). Since dim(50(2, 1)) = 3 (which is the same as the dimension of the 
internal vector space), we can define an 50(2, 1) Lie algebra-valued connection 1-form, 3 A a , 
via 

3 A aI J =: 3 A K a e J IK . (3.27) 

This is just the adjoint representation of the Lie algebra of S0(2, 1) with respect to the 
structure constants e 1 jk '■= f] IM ' £mjk^\ That the defining representation and adjoint repre- 
sentation agree is a property that holds only in 2+1 dimensions since dim(SO(n, 1)) = n + 1 
if and only if n = 2. In terms of 3 A^, the generalized derivative operator ^Da satisfies 

^Dav 1 = d a v I + [ 3 A a ,v] 1 , (3.28) 

where p-A^f] 7 := e 1 jk 3 A J a v K . From (3.27), it also follows that the Lie algebra valued- 
curvature tensor 3 F^ b (which is related to 3 F abI J via 3 F abI J = 3 F^e J jk) can be written as 

3 Fl h = 2d [a 3 Al ] + [ 3 A a , 3 A b ] 1 . (3.29) 

Thus, in terms and 3 Fl b , the Palatini action becomes 

S P ( 3 e, 3 A) = \ [ V abc eiJK % 3 F bc JK 
4 Jm 

= 7 / V abc eiJK V a 3 F b L c e KJ L (3.30) 

4 JM 

f ^.abc 3_ Z-pI 

-2Jm V 6aI Fb - 

But now note that the last line above suggests a natural generalization. Namely, let G 
be any Lie group with Lie algebra £c an d let 3 A T a and 3 e a j be Cq- and £j.-valued 1-forms, 
respectively. Although the action given by (3.30) was originally defined for the Lie group 
5*0(2, 1), it is well-defined in the above sense for any Lie group G. 3 F^ b is still the curvature 
tensor of 3 A ! a , but 3 e a j can no longer be thought of as a co-triad. In fact, since G is now 

10 Given a Lie algebra C with structure constants C 1 jk, the adjoint representation of C by linear op- 
erators on C is defined by the mapping v 1 S C i— > (ad v )i J := v k C j ik- Under ad, the Lie bracket 
[i^w;] 7 := C 1 jkv j w k £ C maps to the commutator of linear operators [ad v , ad w ]i J :— (ad v ) i K (ad w ) k J — 
{ad w ) i K (ad v ) k J ■ I should note that since (ad v )j J w I = — [v,w] J , the above definition of the adjoint rep- 
resentation differs in sign from that given in most math and physics textbooks. The sign difference can be 
traced to my definition of the commutator of linear operators, which also differs in sign from the standard 
definition. 
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arbitrary, the index / can take any value 1, 2, • • • , dim(G). Nonetheless, we can still define 
the Palatini action based on G via 

G S P (\M):=\ [ fj abc3 e a i 3 F^ (3.31) 
2 Jm 

which we treat it as a functional of an /^-valued connection 1-form 3 A ! a and an £^-valued 
covector field 3 e a j. The equations of motion we obtain by varying 3 e aI and 3 A T a are 

~a6c3 F / =0 &nd Sp^-abc 3^ = Q> (3 33) 

which are the analogs of equations (3.16) and (3.17). As before, the second equation re- 
quires a torsion-free extension of to spacetime tensor fields, but again, all results will be 
independent of this choice. 

3.2 Legendre transform 

Given the action (3.31), it is a straightforward exercise to put the 2+1 Palatini theory 
based on G in Hamiltonian form. We will assume that M is topologically £ x R and that 
there exists a function t (with nowhere vanishing gradient (dt) a ) such that each t = const 
surface S t is diffeomorphic to S. As usual, t a will denote the flow vector field satisfying 
t a (dt) a = 1. Since the Lie group G is arbitrary, the 2+1 Palatini theory based on G is not a 
theory of a spacetime metric; it does not involve a spacetime metric in any way whatsoever. 
Thus, in particular, t does not necessarily have the interpretation of time. Nonetheless, we 
can still define "evolution" from one t = const surface to the next using the Lie derivative 
with respect to t a . 

To write (3.31) in 2+1 form, we decompose 

~abc in 

terms of t a and rf h (the Levi-Civita 
tensor density of weig ht +1 on E). Using fj abc = 3t [a rj bc Ut, we get 

b P { e, A) — - n e aI t hc 

Z JM 

= \\ dt + tbf}ca + f?,ab) 3eai ^ (3,33) 
= l dt L \ { " e ■ + ^ A ' c _ ^ Vc{3A ■ t)/? 

where ( 3 e • t)/ := t a 3 e aI , Ef := ij ab 3 e bI , ( 3 A ■ t) 1 := t a 3 Ai, and A* a :— t b a 3 A{ are the 
configuration variables which specify all the information contained in the field variables 3 e a j 
and 3 A{. Note that: 

1. Since G is an arbitrary Lie group, the internal index / can take any value / = 
1,2, • • - ,dim(G). Thus, Ef cannot in general be interpreted as a dyad. In fact, this 
is true even when G = SO(2, 1), since dim(5'0(2, 1)) = 3. However, for SO(2, 1) we 
have EfE bI = \ ab (= qq ab ). 
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2. t a 3 F^ b = C^ 3 A[ — ^Dfy^A ■ t) 1 , which follows from a generalization of Cartan's identity 
Cua = irfda + d(i$a). The Lie derivative treats fields with only internal indices as 
scalars. 

3. Cftl = 0, where t b := 5 b — t a (dt)b is the natural projection operator into the t = const 
surfaces defined by t and t a . 

4. T> a := t a ^Dfe is the generalized derivative operator on E associated with A^. 

5. F* b := t c a t b 3 F^ d is the curvature tensor of V a and satisfies = 2d[ a Al + [A a , A b Y . 

From (3.33), we see that (modulo a surface integral) the Lagrangian G Lp of the 2+1 
Palatini theory based on G is given by 

G L P = jf i( 3 e • t)rrj ab F^ b + EfCf A{ + (V a E^)( 3 A • t) 1 . (3.34) 

By inspection, we see that the momentum conjugate to A a is Ef, while ( 3 e • t)i and ( 3 A ■ t) 1 
both play the role of Lagrange multipliers. Thus, the Dirac constraint analysis says that the 
phase space ( G T P , G Qp) is coordinatized by the pair (A a , Ef) and has symplectic structure^ 

G n P = f &Ef A &Ai. (3.35) 

The Hamiltonian is given by 

G H P (A, E) = jf - l -( 3 e ■ Oi^Fi - (V a Ef)( 3 A ■ t) 1 . (3.36) 

As we shall see in the next subsection, this is just a sum of 1st class constraint functions 
associated with 

fj ab F^ b ^0 and V a E%t*0. (3.37) 

Note that these equations are the field equations (3.32) pulled-back to S with rj ab . Note 
also that they are polynomial in the canonically conjugate variables (A a , Ef). This is to be 
contrasted with the constraint equations for the standard Einstein-Hilbert theory. Recall 
that the scalar constraint of that theory depended non-polynomially on q ab . 

3.3 Constraint algebra 

As usual, to evaluate the Poisson brackets of the constraints and to determine the motions 
they generate on phase space, we must first construct constraint functions associated with 

11 Note that in terms of the Poisson bracket { , } defined by G £lp, we have {A^(x), Ej(y)} — 5 h a 5 I J 5{x : y). 
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(3.37). Given test fields v 1 and aj, which take values in the Lie algebra C G and its dual C* G , 
we define 

F(a) := \ ^ajrj ab F* b and G(v) := jf (3.38) 

We will call G(y ) the Gauss constraint function since it will play the same role as the 
Gauss constraint of Yang-Mills theory. We will see that G(v ) generates the usual gauge 
transformations of the connection 1-form A ! a and it conjugate momentum (or "electric field") 



Ef. 

We are now ready to evaluate the functional derivatives of F(a) and G(v). Since F(a) 
is independent of the momentum Ef, and since 5Fl b = 2V[ a 5Al^ we find 

and ^ = (3.39) 

Similarly, if we vary G{v ) with respect to Ef and A^, we find 

6G{V) = ~Vy and = {v, E«h (:= C^Efi . (3.40) 



SEf 5Ai 

Here C 1 JK denote the structure constants of the Lie algebra C G and { , } : C G x C* G — > C* G 
denotes the co- adjoint bracket. { , } is defined in terms of the Lie bracket [ , ] : C G x C G — > 
C G via {v^^jw 1 := cxk[v,w] k . 

Given (3.39) and (3.40), we can now write down the Hamiltonian vector fields Xp( a ) and 
X G (v) associated with F(a) and G(v). They are 

X F[a) = [ -v a \V bai )-^ and (3.41a) 

Thus, under the 1-parameter family of diffeomorphisms on G Tp associated with X F ^, we 
have 

Ai ^ Ai and (3.42a) 
Ef i-> Ef -eirfV^) +0(e 2 ). (3.426) 

Similarly, under the 1-parameter family of diffeomorphisms on G T P associated with X G ^, 
we have 

A l~ + 0(e 2 ) and (3.43a) 

Ef ^Ef- e{v, E a }! + 0(e 2 ). (3.436) 
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Note that (3.43a) and (3.436) are the usual gauge transformations of the connection 1-form 
and its conjugate momentum Ef that we find in Yang-Mills theory. Equations (3.42) and 
(3.43) are the motions on phase space generated by F(a) and G(v). 

Given (3.39) and (3.40), we can also evaluate the Poisson brackets of the constraint 
functions. Since the G(v )'s play the same role as the Gauss constraint functions of Yang- 
Mills theory, we would expect their Poisson bracket algebra to be isomorphic to the Lie 
algebra Cc- This is indeed the case. We find 

{G(v),G(w)} = G([v,w]), (3.44) 

where [f,w] 7 = C 1 jkv j w k is the Lie bracket of v 1 and w 1 . Thus, v 1 i— > G(v) is a represen- 
tation of the Lie algebra Cg- The Lie bracket in Cg is mapped to the Poisson bracket of the 
corresponding Gauss constraint functions. 

Since F(a) is independent of Ef, it follows trivially that 

{F(a),F(f3)} = 0. (3.45) 

With only slightly more effort, we can show that 

{G(v),F(a)} = -F({v,a}), (3.46) 

where {v, a}/ = C K jiv J ax is the co-adjoint bracket of v 1 and a/. Thus, the totality 
of constraint functions (G(v) and F(a)) is closed under Poisson bracket — i.e., they form 
a 1st class set. Furthermore, since (3.44), (3.45), and (3.46) do not involve any structure 
functions (unlike the constraint algebra of the Einstein-Hilbert theory discussed in subsection 
2.3), the set of constraint functions form a Lie algebra with respect to Poisson bracket. 
In fact, (a, v) G C* G x C G \— > (F(a),G(v)) is a representation of the Lie algebra Cjg of 
the inhomogeneous Lie group IG associated with G.Q The action T v {a) := — {v,a}i of 
v 1 G Cg on aj G C* G is mapped to the Poisson bracket {G(v), F(a)} of the corresponding 
constraint functions. The F(cn)'s play the role of "translations" and the G(v)'s play the role 
of "rotations" in the inhomogeneous group. 

4. Chern-Simons theory 

So far in this review, we have written down two different actions for 2+1 gravity: The 
standard Einstein-Hilbert action, which gave us a description of 2+1 gravity in terms of a 

12 We will discuss the construction of the inhomogeneous Lie group IG and its Lie algebra Cig in subsection 
4.4 where we show the equivalence between the 2+1 Palatini theory based on G and Chern-Simons theory 
based on IG. When G is the Lorentz group SO(2, 1), IG is the corresponding Poincare group ISO(2, 1). 
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spacetime metric (or equivalently, a co-triad); and the 2+1 Palatini action based on £0(2, 1), 
which gave us a description in terms of a co-triad and a connection 1-form. In this section, 
we shall see that 2+1 gravity can be described by an action that depends only on a con- 
nection 1-form. We shall see that the 2+1 Palatini action based on £0(2, 1) is equal to the 
Chern-Simons action based on IS0(2, 1), modulo a surface integral that does not affect the 
equations of motion. This result was first shown by A. Achucarro and P.K. Townsend 
it was later rediscovered and used by Witten || to quantize 2+1 gravity. In this section, 
we will follow the treatment of |SD| in which the result for 2+1 gravity follows as a special 
case. We will show that the 2+1 Palatini theory based on any Lie group G is equivalent to 
Chern-Simons theory based on the inhomogeneous Lie group IG associated with G. 

The above equivalence between the 2+1 Palatini and Chern-Simons theories is at the 
level of actions. The gauge groups, G and IG, are different, but the actions are the same. 
It is interesting to note that the 2+1 Palatini and Chern-Simons theories based on the same 
Lie group G are also related, but this time at the level of their Hamiltonian formulations. We 
shall see that the reduced phase space of the Chern-Simons theory based on G is the reduced 
configuration space of the 2+1 Palatini theory based on the same G. Since Chern-Simons 
theory is not available in 3+1 dimensions, the relationships that we find in this section do 
not, unfortunately, extend to 3+1 theories of gravity. 



4.1 Euler-Lagrange equations of motion 

Unlike the standard Einstein-Hilbert and Palatini theories which are well-defined in n+1 
dimensions, Chern-Simons theory is defined only in odd dimensions. In 2+1 dimensions, the 
basic variable is a connection 1-form 3 A\ which takes values in a Lie algebra Lq equipped 
with an invariant, non-degenerate bilinear form The Chern-Simons action based on G 

is defined by 

G Scs( 3 A) := \ J M f, abc k l3 ( *A% 3 A{ + ~ 3 A\[ 3 A b 3 A c }^ , (4.1) 

where [ 3 A b , 3 A c ] j := C j mn 3 A\ n 3 A n c denotes the Lie bracket of 3 A\ and 3 A l c . It is important 
to note that Chern-Simons theory is not defined for arbitrary Lie groups — we need the 
additional structure provided by the invariant, non-degenerate bilinear form fey. 

To obtain the Euler-Lagrange equations of motion, we vary G Scs( 3 A) with respect to 3 A l a . 

13 We will require that fey be invariant under the adjoint action of the Lie algebra Lq on itself — i.e., that 
kij[x,v] l wi + kijV l [x, u>p = for all v % ,w l ,x l € Cq- If Cijk is defined in terms of the structure constants 
C l jk via Cijk '■= ki m C m jk, then invariance of fcy under the adjoint action is equivalent to Cijk = C[yfc]- If 
the Lie group is semi-simple (i.e., if it does not admit any non-trivial abelian normal subgroup), then we 
are guaranteed that such a fey exists. This is just the Cartan- Killing metric defined by k%j := C m n iC n m j- 
Invariance of is equivalent to the invariance of C l jk — that is, the Jacobi identity C m [ijC n w m = 0. 
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Using the fact that := k im C m jk is totally anti-symmetric, we obtain 

= 0, (4.2) 

where 3 F* b = 2d[ a 3 A^ + [ 3 A a , 3 A b Y is the Lie algebra- valued curvature tensor of the generalized 
derivative operator ^Da defined by 

^■^darf + i^vy. (4.3) 

If we also use the fact that fc^ is non-degenerate, we get 3 F^ b = 0. Thus, Chern-Simons 
theory is a theory of a flat connection 1-form. We will see the role that this equation plays 
in the next two subsections when we put the theory in Hamiltonian form. 

4.2 Legendre transform 

Just like the 2+1 Palatini theory based on an arbitrary Lie group G, Chern-Simons theory 
is not a theory of a spacetime metric. However, we can still put this theory in Hamiltonian 
form if we assume that M is topologically E x R for some sub manifold E and assume that 
there exists a function t (with nowhere vanishing gradient (dt) a ) such that each t = const 
surface E t is diffeomorphic to E. As usual, we let t a denote the flow vector field satisfying 
t a (dt) a = 1. 

Given t and t a , we are now ready to write the Chern-Simons action (4.1) in 2+1 form. 
Using the decomposition rj abc = 3t^ a fj bc ^dt, we get 

G S CS ( 3 A) = \f *f%( 3 A% 3 Ai + \ 3 A\[ 3 A b , 3 A c f) 
= \jdt jj*A ■ t)%rf^Fi + fj^hjAiLr Ai, 

where the last equality holds modulo a surface integral. Here ( 3 A ■ t) 1 := t a 3 A l a and A\ := 
t b a 3 A\ (= (5 b a — t b (dt) a ) 3 A\) are the configuration variables which specify all the information 
contained in the field variable 3 A l a . The Lie derivative treats fields with only internal indices 
as scalars, and F % ah = 2d[ a A^ + [A a , A b } 1 is the curvature tensor associated with A\. 

From (4.4), it follows that the Lagrangian G L C s of the Chern-Simons theory based on G 
is given by 

G Lcs = I J/A ■ t)%tf b Fi h + ^MiA- K- (4-5) 

The momentum canonically conjugate to A l a is ^fj ab kijA J b , while ( 3 A ■ t) z plays the role 
of a Lagrange multiplier. Thus, the Dirac constraint analysis says that the phase space 
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(4-4) 



( g Tqs, g Qcs) is coordinatized by (A l a ) and has symplectic structure^ 

1 



G 



The Hamiltonian is given by 



G H CS (A) = -\ J ('A ■ tYk^Fl (4.7) 



As we shall see in the next subsection, this is just a 1st class constraint function associated 
with 

hiV ab Fi b = 0. (4.8) 

Note that constraint equation (4.8) is the field equation (4.2) pulled-back to S with rj ab . Just 
as in the 2+1 Palatini theory, the constraint equation is polynomial in the basic variable A\. 
Note also that although (4.8) may not look like the standard Gauss constraint of Yang-Mills 
theory, we shall see that its associated constraint function generates the same motion of A l a . 

4.3 Constraint algebra 

Following the same procedure that we used in Sections 2 and 3, we first construct a 
constraint function associated with (4.8). Given a test field v % (which takes values in the Lie 
algebra Co), we define 

G{v) :=)-( v^ b Fi b . (4.9) 



2 7s 

Since the phase space is coordinatized by the single field A l a , we need to evaluate only one 
functional derivative. Varying G(v ) with respect to A l a , we get 

,)(T ' (r) -k^ ab W, (4.10) 



5A 



a 



where T> a is any torsion-free extension of the generalized derivative operator associated with 
A\. From (4.10) it then follows that the Hamiltonian vector field Xq{ v ) is given by 

X G(v)= L-( V aV l )j^, (4.11) 



a 



so that 

A\^A\-eV a v* + 0{e 2 ) (4.12) 



14 Note that in terms of the Poisson bracket { , } defined by G f2cs, we have {A % a (x) , A J b (y)} = r^bk^ 5(x, y), 
where r^b and fc u denote the inverses of rj ab and kij. This result follows from the fact that for any / : G Tcs — > 

/ Sg 



R, the Hamiltonian vector field Xf is given by Xf = jj.flpbk^ jjj jjr- Hence the Poisson bracket of any 
two functions /, g : G r cs - R is {/, g} = / E j^k** fa jh- 
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under the 1-parameter family of diffeomorphisms on Tcs associated with Xqm- This is 
the usual gauge transformation of the connection 1-form that we find in Yang-Mills theory. 
Thus, G(v) can be appropriately called a Gauss constraint function. 

Given (4.10), it is also straightforward to evaluate the Poisson brackets of the constraints. 
We find that 

{G(v),G(w)} = G([v 1 w}), (4.13) 

which is the expected Poisson bracket algebra of the Gauss constraint functions. The map 
v l \— > G(v) is a representation of the Lie algebra Cq- The Lie bracket in C G is mapped to 
the Poisson bracket of the corresponding Gauss constraint functions. 

4.4 Relationship to the 2+1 Palatini theory 

Before we can show the relationship between the Chern-Simons and 2+1 Palatini theories, 
we will first have to recall the construction of the inhomogeneous Lie group IG associated 
with any Lie group G. This will allow us to generalize the equivalence of the 2+1 Palatini 
and Chern-Simons theories (as shown in || |32f ) to arbitrary Lie groups G. We will be able 
to show that the 2+1 Palatini theory based on any G is equivalent to Chern-Simons theory 
based on IG. 

Consider any Lie group G with Lie algebra Cg, and let C* G denote the vector space dual 
of Cg- If v 1 , w 1 denote typical elements of Cq and ckj, (3j denote typical elements of C G , 
then (a, v) 1 := (a/,f / ) and (/3,w) 1 := (/3j, w 1 ) are typical elements of the direct sum vector 
space C* G © Cq- We can define a bracket on C* G © Cg via 

[(a, v), (0, w)}' := (-{v, (3} + {w, a}, [v, w})\ (4.14) 

where [f,w] 7 := C i jkv j w k and {v,/3}j := C K jiv j (3k are the Lie bracket and co-adjoint 
bracket associated with Lq- By inspection, we see that (4.14) is linear and anti-symmetric. 
If we use 

{[v, w],a}j = -{v, {w, a}}/ + {w, {v, (4.15) 

(which follows as a consequence of the Jacobi identity for Cg), we can show that (4.14) 
satisfies the Jacobi identity as well. Thus, the vector space Cjg '■= £-h ® ^-g together with 
(4.14) is actually a Lie algebra. We call Cjg the inhomogeneous Lie algebra associated with 
G; the inhomogeneous Lie group IG is obtained by exponentiating the Lie algebra Cjg- As 
we shall see later in this subsection, IG is simply the cotangent bundle over G. 

The terminology inhomogeneous is due to the fact that Cjg admits an abelian Lie ideal 
isomorphic to C Gl and that the quotient of Cig by this ideal is isomorphic to £g-0 Thus, 

15 A Lie ideal I of a Lie algebra £ is a vector subspace X C C such that [i, i] el for any i £l, x G C. 
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elements of C* G are analogous to infinitesimal "translations," while elements of Lq are analo- 
gous to infinitesimal "rotations." Note, however, that the space of translations and rotations 
have the same dimension. As a special case, if one chooses G to be the 2+1 dimensional 
Lorentz group 5*0(2, 1), then the above construction yields for IG the 2+1 dimensional 
Poincare group 750(2, 1). 

In addition to the above Lie algebra structure, C* G © Cq is equipped with a (natural) 
invariant, non-degenerate bilinear form kij defined by 



Since Lie is not semi-simple (because it admits a non-trivial abelian Lie ideal), is not the 
(degenerate) Cartan-Killing metric of Ciq. Nevertheless, the existence of will allow us 
to construct Chern-Simons theory for IG. Recall that without an invariant, non-degenerate 
bilinear form, Chern-Simons theory could not be defined. Note also that for G = 5*0(2, 1), 
the above construction of k^ reduces to that used by Witten ||. 

Given these remarks, we can now show that the 2+1 Palatini theory based on any Lie 
group G is equivalent to Chern-Simons theory based on IG. To do this, recall that for any 
Lie group G, the 2+1 Palatini action based on G is given by 



where 3 A^ and 3 e a i are Cq- and ^-valued 1-forms. We now construct the inhomogeneous 
Lie algebra Cjg associated with G and define an £/G-valued connection 1-form 3 A l a via 



By simply substituting this expression for 3 A l a into the Chern-Simons action IG S C s( 3 A), we 



kij{a,v)\(3,wy := ajw 1 + ftiv 1 . 



(4.16) 




(4.17) 



3 A l a '■— ( 3 e a /, 3 -^a)- 



(4.18) 
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find that 

IG 



(4.19) 



; SgsCA) = \ J M V abc ki( 3 A\d h 3 A{ + \ "A\['A h *A c \ j ) 
= \l V ab \"e aI d b ^ c + {d b \ cI fA{ 

1 JM x 

+ \{ 3 e aI [ 3 A b , *A C Y - {%, 3 e c }, 3 ^ + {*A C , %}, 3 ^)) 
= \l V abc ( 3 e a id b 3 Al + d^eJK) - 3 e cI d b 3 A{ 

+ \( 3 e a i[ 3 A b , 3 A C Y - 3 e c7 [X 3 A a Y + \ bI ?A c , 3 A a Y)) 
= \J m V abc { 3 e a i(2d b 3 Ai) + 3 e aI [ 3 A b , *A C ]> + d b ( 3 e c Mi)) 
= - [ r\ ahc 3 e aI 3 F/ C + (surface integral) 

2 JM 

= G S P ( 3 e, 3 A) + (surface integral), 

where we have used definitions (4.14), (4.16), and (4.18) repeatedly. Since the surface term 
does not affect the Euler-Lagrange equations of motion, we can conclude that the 2+1 
Palatini theory based on G is equivalent to Chern-Simons theory based on IG. This is the 
desired result. Note that as a special case, we can conclude that 2+1 gravity as described 
by the 2+1 Palatini action based on 5*0(2, 1) is equivalent to Chern-Simons theory based 
on 750(2,1). 

Up to now, we have only described the inhomogeneous Lie group IG in terms of its asso- 
ciated Lie algebra Cjq- We just exponentiated the Lie algebra Cig to obtain IG. However, 
it is also instructive to give an explicit construction of IG at the level of groups and mani- 
folds. But to do this, we will need to make another short digression, this time on semi-direct 
products and semi-direct sums. Readers already familiar with these definitions may skip to 
the paragraph immediately following equations (4.24). 

Let G and H be any two groups. To define the semi-direct product H @ G, we need a 
homomorphism a from the group G into the group of automorphisms of H — i.e., for each 
g,g' G G and h', h' G H, the map a g : H — * H must be 1-1, onto, and satisfy 

a g {hh')=a g {h)a g {h') and a gg ,(h) = <J g {a g ,{h)) (4.20) 

for every g,g' G G and h', hi G H. Given this structure, one can check that 

(h,g)(h',g') = (ha g (h'),gg') (4.21) 

defines a group multiplication law on the set H x G. The identity element is (en, £g), where 
en, £g are the identities in H and G, and the inverse (h.g)' 1 of (h,g) is (o- g -i{h~ v ), g~ v ). 
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The set H x G together with this group multiplication law defines the semi- direct product 
H @ G. Note that H is homomorphic to a (not necessarily abelian) normal subgroup of 
H @ G, and the quotient of H @ G by this normal subgroup is homomorphic to G. As a 
trivial example, if a g {h) = h for all g G G, h G H, then H @ G is the usual direct product 
H eg) G of groups. 

Now assume that G and if are Lie groups with Lie algebras Cg and Ch- If H ©G is the 
semi-direct product of G and if with respect to some action a of G on if satisfying (4.20), 
we would now like to determine the relationship between the Lie algebras Chqg, £>g> an d 
Ch- To do this, we differentiate the action o of G on H to obtain an action r of Cg on 
Ch- More precisely, if g(e) is a 1-parameter curve in G with g(Q) = cg and tangent vector 
v := j-| e=0 <7(e) and h(X) is a 1-parameter curve in H with /i(0) = e# and tangent vector 
a := 4r\ x=0 h(\), then we define 



d 



de 



d 



-o d\ 



-rvipL) ■= T. _ "7T ^)(MA)). (4.22) 



A=0 



In terms of r, the Lie bracket of (a,v) and (fl,w) in Chqg becomes^ 

[(a,u), (/3,w)] = (r„(/3) - r w (a) + [«,«;]), (4.23) 

where [u, w] and [a, 0\ are the Lie brackets of v,w G Co and a,/? 6 Note that (4.23) 
satisfies the Jacobi identity as a consequence of 

r v ([a,p]) = [t v (o),0\ + [a,r v (P)} and (4.24a) 
T[v,w]( a ) = t v{t w {u)) - T w (T v (a)), (4.246) 

which follow from the definition of r and the properties (4.20) satisfied by a. Thus, if Cg 
and Ch are two Lie algebras and r is an action of Cq on Ch satisfying (4.24), then the direct 
sum vector space Ch © Cg together with the bracket defined by (4.23) is a Lie algebra. This 
Lie algebra, denoted Ch ©Cg, is called the semi-direct sum of Cg and Ch- Note that Ch is 
isomorphic to a (not necessarily abelian) Lie ideal of Ch ©Cg, and the quotient of Ch ©Cg 
by this ideal is isomorphic to Cg- If cr g (h) = h for all g G G, h G H (so that H©G = H®G), 
then Ch ©Cg is the usual direct sum Ch © Cg of Lie algebras. 

Given these general remarks, let us now return to our discussion of the inhomogeneous 
Lie group IG and its associated Lie algebra Cjg- From the above definitions, we see that 

16 To obtain this result, consider 1-parameter curves (h(e),g(e)) and (h'(e'),g'(e')) in H @ G with 
(ft(0),s(0)) = (h'(0),g'(0)) = (e H ,e G ) and tangent vectors (a,v) := g~|e=o(M e )> s( e )) and w ) : = 
xr| e '=o(/i / (e')) ff'( e '))- Then use the definition of the Lie bracket in terms of the group multiplication law 



(4.21), [(a.t)),^^)] := £U =0 £U>=o(h(e), g{e)){h'{e'),g'{e')){h{e),g{e))-^h'{e'),g'{e'))-\ This leads to 



(4.23). 
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Cig is simply the semi-direct sum C* G (t)Cg- C* g is to be thought of as a Lie algebra with 
the trivial Lie bracket [a, [3] = for all ai, Pi G C G ; the action r of Co on C* G is given 
by t v (P) = —{v,P}j. Equations (4.24) hold for this action as a consequence of the Jacobi 
identity in C G : Equation (4.24a) is satisfied since [a,P] = for all aj, Pi G C G , while 
equation (4.246) is equivalent to equation (4.15). Furthermore, the inhomogenized Lie group 
IG is simply the semi-direct product C* G @ G. C* G is to be thought of as an abelian group 
with respect to vector addition, and the action a of G on C* G is induced by the adjoint action 
of G on itself.0 This implies that as a manifold IG is the cotangent bundle T*G. At each 
point g G G, the cotangent space T*G is isomorphic to C* G . 

Moreover, the above relationship between G and IG allows us to prove an interesting 
mathematical result involving the space of connection 1-forms on a 2-dimensional manifold. 
We can show that for any Lie group G 

T*( G A) = IG A, (4.25) 

where G A and IG A denote the space of Cg- and £/G-valued connection 1-forms on a 2- 
dimensional manifold S. The map 

(Ai, E}) G T*( G A) ^ A\ := (e aI , A 1 .) G IG A (4.26) 

(where e a j := —r^^Ej) is a diffeomorphism from the manifold T*( G A) to the manifold IG A 
that sends the natural symplectic structure 

G Vl := J &E? A dAi (4.27) 

on T*( G A) to the natural symplectic structure 

IG Q := -- f ?f 6 %dL4 l a A&Ai (4.28) 

on IG A. (Here kij denotes the (natural) invariant, non-degenerate bilinear form on Cig 
defined by (4.16).) Note that (4.27) and (4.28) are the symplectic structures of the 2+1 
Palatini theory based on G and the Chern-Simons theory based on IG. However, the above 
result (4.25) does not require any knowledge of the 2+1 Palatini or Chern-Simons actions. 

Finally, to conclude this subsection, I would like to verify the claim made at the start of 
Section 4 that the reduced phase space of Chern-Simons theory based on any Lie group G 

17 The adjoint action of G on itself is denned by A g (g') — gg'g^ 1 for all g,g' £ G. By differentiating 
A g at the identity e, we obtain a map Ad g : Cg — > Cg via Ad g {v) :— A' g (e) ■ v. Ad defines the adjoint 
representation of the Lie group G by linear operators on the Lie algebra Cg- The action a of G on C* G is 
then given by (o~ g (a))(v) := a(Ad g (v)) for any a € C G and v £ Cg- 
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is the reduced configuration space of the 2+1 Palatini theory based on the same Lie group 
G. This result will be simpler to prove than the previous two results since most of the 
preliminary work has already been done. 

Let G be any Lie group whose Lie algebra Lq admits an invariant, non-degenerate bilinear 
form kjj. Then Chern-Simons theory based on G is well-defined, and, as we saw in subsection 
4.2, the phase space G Tcs is coordinatized by /^-valued connection 1-forms A ! a on S. The 
symplectic structure is 

G n cs = ~ I rj^ku&Ai A dL4 6 J . (4.29) 
2 js 

In subsection 4.3, we then verified that the constraint functions G(v) associated with the 
the constraint equation 

kjj^Fi = (4.30) 
formed a 1st class set and generated the usual gauge transformations 

Ai^Al-eVy + Oie 2 ). (4.31) 

To pass to the reduced phase space, we must factor-out the constraint surface (defined by 
(4.30)) by the orbits of the Hamiltonian vector fields X^)^ From (4.30) and (4.31) we see 
that the reduced phase space G Tcs of the Chern-Simons theory based on G is coordinatized 
by equivalence classes of flat Lq- valued connection 1-forms on S, where two such connection 
1-forms are said to be equivalent if and only if they are related by (4.31). This space is called 
the moduli space of flat Lq- valued connection 1-forms on E. 

Now recall the Hamiltonian formulation of the 2+1 Palatini theory based on the same 
Lie group G. In subsection 3.3, we saw that the phase space G Tp was coordinatized by pairs 
(Aq, Ej) consisting of /^-valued connection 1-forms A\ on £ and their canonically conjugate 
momentum Ef. G Tp was the cotangent bundle T*( G Cp) over the configuration space G Cp of 
ZZc-valued connection 1-forms A\ on E with symplectic structure 

G n P = f &Er A dAi. (4.32) 

We also saw that the constraint equations of the 2+1 Palatini theory were 

V ab F^ b = and V a Ef = 0, (4.33) 

18 Recall that given a symplectic manifold (r, fi), a set of constraints <fii form a 1st class set if and only if 
each Hamiltonian vector field X ( p i is tangent to the constraint surface V C V defined by the vanishing of all 
the constraints. The pull-back, il, of fl to T is degenerate with the degenerate directions given by the X r j }i . 
Thus, (r, f2) is not a symplectic manifold. However, by factoring-out the constraint surface by the orbits of 
the Xfa, we obtain a reduced phase space (r, fl) whose coordinates are precisely the true degrees of freedom 
of the theory. Cl is non-degenerate; it is the projection of fl to T. 
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and verified that their associated constraint functions F(a) and G(v) formed a 1st class set. 
They generated the motions 



Ai^Ai and (4.34a) 

Ef i-> El - eirfV^j) + 0(e 2 ) (4.346) 

and 

A{ >-> A{ - eVaV 1 + 0(e 2 ) and (4.35a) 

Ef ^ -e{v,£ a } 7 + 0(e 2 ), (4.356) 

respectively. Thus, the reduced phase space G Tp of the 2+1 Palatini theory based on G is 
coordinatized by equivalence classes of pairs consisting of flat £c-valued connection 1-forms 
on E and divergence- free C* G - valued vector densities of weight +1 on E, where two such pairs 
are said to be equivalent if and only if they are related by (4.34) and (4.35). Since F(a) 
and G(v) are independent and linear in the momentum Ef, it follows that G Tp is naturally 
the cotangent bundle T*{ G Cp) over the reduced configuration space G Cp of the 2+1 Palatini 
theory based on G. From (4.34a) and (4.35a) we see that G Cp is again the moduli space of 
flat /^-valued connection 1-forms on E. Thus, G t C s — G Cp as desired. In particular, the 
reduced configuration space of the 2+1 Palatini theory based on G has the structure of a 
symplectic manifold. 

This last result has interesting consequences. It can be used, for example, to show the 
relationship between the T°[j] and T 1 ^] observables for 2+1 gravity. These are the 2+1 
dimensional analogs of the classical T-observables constructed by Rovelli and Smolin for 
the 3+1 theory. As shown in |3(J, T°[j] is the trace of the holonomy of the connection 



around a closed loop 7 in E, while ^[7] is the function on the reduced phase space of the 
2+1 Palatini theory defined by the Hamiltonian vector field associated with T°[7]. Thus, 
many properties satisfied by the T^J's can be derived from similar properties satisfied by 
the T°[7]'s. 



5. 2+1 matter couplings 

In this section, we will couple various matter fields to 2+1 gravity via the 2+1 Palatini 
action. We will consider the inclusion of a cosmological constant and a massless scalar field. 
One can couple other fundamental matter fields (e.g., Yang- Mills and Dirac fields) to 2+1 
gravity in a similar fashion — I have chosen to consider a massless scalar field in detail since 
2+1 gravity coupled to a massless scalar field is the dimensional reduction of 3+1 vacuum 
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general relativity with a spacelike, hypersurface-orthogonal Killing vector field [16|. As noted 
in Section 1, this is an interesting case since it appears likely that the non-perturbative 
canonical quantization program for 3+1 gravity can be carried through to completion for 
this reduced theory. 

In subsection 5.1, we define the 2+1 Palatini theory based on a Lie group G with cos- 
mological constant A. We derive the Euler-Lagrange equations of motion and perform a 
Legendre transform to obtain a Hamiltonian formulation of the theory. Just as in Section 
3 (when A was equal to zero), we shall find that the constraint equations of the theory are 
polynomial in the canonically conjugate variables. We shall also find that their associated 
constraint functions still form a Lie algebra with respect to Poisson bracket. 

In subsection 5.2, we will show that the 2+1 Palatini theory based on G with cosmological 
constant A is equivalent to Chern-Simons theory based on the A-deformation, AG, of G. This 
is a generalization of Witten's result § for G = SO(2, 1) (and AG = SO(3, 1) or SO(2, 2) 
depending on the sign of A) which holds for any Lie group G that admits an invariant, totally 
anti-symmetric tensor e IJK . This result also generalizes the A = equivalence of the 2+1 
Palatini and Chern-Simons theories given in subsection 4.4. 

Finally, in subsection 5.3, we define the action for a massless scalar field and couple this 
field to 2+1 gravity by adding the action to the 2+1 Palatini action based on G = SO(2, 1). 
It is when we wish to couple matter with local degrees of freedom to 2+1 gravity (as it is 
for the case of a massless scalar field) that we are forced to take the Lie group G to be 
such that the fields 3 e^ have the interpretation of a co-triad. We obtain the Euler-Lagrange 
equations of motion for the coupled theory and then perform a Legendre transform to obtain 
the Hamiltonian formulation. We shall find that the constraint equations remain polynomial 
in the canonically conjugate variables and the associated constraint functions form a 1st 
class set, but they no longer form a Lie algebra with respect to Poisson bracket. 

The basis for much of the material in this section can be found in || |30], ^ . 



5.1 2+1 Palatini theory coupled to a cosmological constant 

Recall that the equation of motion for gravity coupled to the cosmological constant A is 

G ab + Ag ab = 0, (5.1) 

where G a b '■= Rab — \Rg a b is the Einstein tensor of g a b- We can obtain (5.1) via an action 
principle if we modify the standard Einstein-Hilbert action by a term proportional to the 
volume of the spacetime. Defining 



SA(g ab ):= I 2A), (5.2) 
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we find that the variation of (5.2) with respect to g ab yields (5.1) (modulo the usual boundary 
term associated with the standard Einstein-Hilbert action). These results are valid in n+1 
dimensions. 

To write this action in 2+1 Palatini form, we proceed as in Section 3. We replace 
the spacetime metric g ab with a co-triad 3 e a i and replace the unique, torsion-free spacetime 
derivative operator V a (compatible with g ab ) with an arbitrary generalized derivative operator 
'*D a . Recalling that ■ s f zz g~ = ^rj abc e IJK 3 e a j 3 e b j 3 e C K, we define 

Q /'3_ 3a\ . _ 1 f ^abc 3_ -irpl ^ f -abcJJK 3_ Z c 3„ 

oa{ e, A) . — - r] e aI b bc - — / r] e e aI e bJ e cK 

2 Jm 3! Jm oX 

(5.3) 



1 f ~abc 3_ / 3 P / ^ JJK 3_ 3_ \ 

- / V e ai{ -rb c -^-e e bJ e cK ), 

2 Jm 3 



where 3 Fl b = 2d[ a 3 A^ + [ 3 A a , 3 A b Y is the internal curvature tensor of the generalized derivative 
operator defined by 

^t/ :=dy +[ 3 A a ,vY. (5.4) 

Note that [ 3 A a , v} 1 := e 1 jk 3 A J a v K where e 1 jk '■= ^ imn Vmj t Ink- Just as we did for the vacuum 
2+1 Palatini theory, we have included an additional overall factor of 1/2 in definition (5.3). 

Although the action (5.3) was originally defined for G = 5*0(2, 1), it is well-defined for 
any Lie group G that admits an invariant, totally anti-symmetric tensor e /Ji<: .|^| This is 
additional structure that does not naturally exist for an arbitrary Lie group G, so unlike 
the 2+1 Palatini theory with A = 0, the 2+1 Palatini theory with non-zero cosmological 
constant A is not defined for arbitrary G. If the Lie algebra Cq admits an invariant, non- 
degenerate bilinear form kjj, then we are guaranteed that such an e IJK exists — we can 
take e IJK := k JM k KN C I mn- Thus, in particular, 2+1 Palatini theory with a non-zero 
cosmological constant is well-defined for semi-simple Lie groups. We should emphasize, 
however, that it is not necessary to restrict ourselves to semi-simple Lie groups. In what 
follows, we will only assume that e IJK exists. Given such a Lie group G, the action 

%( 3 e, 3 A) := l -( f} abc 3 e aI ( 3 F b { - ±e IJK 3 e bJ 3 e cK ) (5.5) 
I JM A 

will be called the 2+1 Palatini action based on G with cosmological constant A. Note that 

3 F I [ b = 2d [a 3 A I b] + i 3 A a , 3 Ab} 1 , (5.6) 

19 We will require that e IJK be invariant under the adjoint action of the Lie algebra Cq on its dual L* G — 
i.e., that e IJK {v, a}i/3jjK + e IJK ai{v, 0}jjk + e IJK ai(3j{v 7 "f}K = for all v 1 € Cg and a/, /3j, 77 G 
C G . ({v,a}i is the co-adjoint bracket of v 1 and ai defined in terms of the structure constants C 1 jk via 
{v,a}i :— C K jiV j olk-) Invariance of e IJK is equivalent to e M ^ IJ C k \in = 0. 
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where [ 3 A a , 3 A b ] ! := C 1 JK 3 A J a 3 Af is the Lie bracket in C G . It is only for G = SO(2, 1) that 

^ JK — £ JK — £ VMJVNK- 

To obtain the Euler-Lagrange equations of motion, we vary G S < A( 3 e, 3 A) with respect to 
both 3 e aI and 3 ,4^. We find 

~ ab c { 3p£ _ Ae IJK 3 6ck) = and Sp^c 3^ = 0; (5 J) 

where the second equation, as usual, requires a torsion-free extension of the generalized 
derivative operator to spacetime tensor fields, but is independent of this choice. Note 
further that if A 7^ 0, is not flat. In fact, for the special case G = SO(2, 1), equations 
(5.7) imply that the spacetime (M,g ab := 3 e a i 3 e b ji] IJ ) has constant curvature equal to 6A. 

To show that the above two equations reproduce (5.1), let us restrict ourselves to G — 
SO(2, 1) (with e IJK being the volume element of rju) so that 3 e aI is, in fact, a co-triad. 
Then by following the argument given in Section 3, we find that the second equation implies 
3 A T a = T T a , where T T a is the (internal) Christoffel symbol of the unique, torsion-free generalized 
derivative operator V a compatible with the co-triad 3 e a j. Thus, is not arbitrary, but 
equals V a when acting on internal indices. Substituting this solution back into the first 
equation, we find 

V abc (R I bc -Ae IJK3 e bJ 3 e cK ) = 0, (5.8) 
where R 1 ^ is the (internal) curvature tensor of V a . Contracting (5.8) with 3 ej gives 

Q ad + Ag ad = Q ( 59 j 

This is the desired result. 

To put this theory in Hamiltonian form, we will assume that M is topologically E x R, 
and assume that there exists a function t (with nowhere vanishing gradient (dt) a ) such that 
each t = const surface £4 is diffeomorphic to £. By t a we will denote the flow vector field 
satisfying t a (dt) a = 1. Using rj abc = 3t^ a rj bc ^dt (and our decomposition of G Sp( 3 e, 3 A) from 
Section 3), we obtain 

G S A ( 3 e, 3 A) = fdtf \{ 3 e ■ t)r(rj^ b - Ae'^E^) 

j z (5.10) 

+ E a I C T A I a -E a I V a { 3 A-t) 1 . 

The configuration variables are ( 3 e • t)j := t a 3 e a i, Ef := rj ab 3 e b i, ( 3 A ■ t) 1 := t a 3 A £ a , and 
A ! a := t b a 3 A[. Thus, (modulo a surface integral) the Lagrangian G L\ of the 2+1 Palatini 
theory based on G with cosmological constant A is given by 

G Lk = [ k 3 e ■ t) T (r b FL - Ae^E^) 

JS z (5.11) 

+ EfC t Ai + (V a Ef)( 3 A.t) 1 . 
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g La is to be viewed as a functional of the configuration variables and their first derivatives. 

Following the standard Dirac constraint analysis, we find that the momentum canonically 
conjugate to A T a is Ef, while ( 3 e • t)j and ( 3 A ■ t) 1 both play the role of Lagrange multipliers. 
Thus, the phase space and symplectic structure are the same as those found for the 2+1 
Palatini theory with A = 0, and the Hamiltonian is given by 

G H A (A, E) = jf -i( 3 e • tUfj^F^ - Ae IJK r^ b E a jE b K ) - (V a Ef)( 3 A ■ t) 1 . (5.12) 

We will see that this is just a sum of 1st class constraint functions associated with 

rj ab F* b - Ae IJK r^ b E a jE b K «0 and D a £J«0. (5.13) 

By inspection, constraint equations (5.13) are polynomial in the canonically conjugate vari- 
ables (A^Ef). They are the field equations (5.7) pulled-back to £ with rj ab . 

As usual, given test fields ai and v 1 , which take values in C* G and Cq, we can define 
constraint functions 

F(a) := \ jf a^Fl, - Ke IJK r^ b E a jE b K ) and G{v) := ^ W 7 (I> a S?). (5.14) 

Note that G(v) is unchanged from the 2+1 Palatini theory with A = 0, while F(a) has an 
additional term quadratic in the momentum Ef. There is only one new functional derivative, 

= -Ae^ b E b a K . (5.15) 

All the others are the same as before. 

Under the 1-parameter family of diffeomorphisms associated with the Hamiltonian vector 
field Xp( a ), we have 

A T a h^ A* a — e(Ae IJK r] slb Eja K ) + 0(e 2 ) and (5.16a) 
E a j^E a j- e(rj ab V bai ) + 0(e 2 ). (5.166) 

Similarly, under the 1-parameter family of diffeomorphisms associated with the Hamiltonian 
vector field Xq( v ), we have 

A ! a h-> A T a - eVy + 0(e 2 ) and (5.17a) 
E?t-+ Ej -e{v,E a }i + 0(e 2 ). (5.176) 

Comparing these results with those from the 2+1 Palatini theory with A = 0, we see that the 
motion of 3 A^ generated by the constraint functions no longer corresponds to the usual gauge 
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transformation of Yang-Mills theory. This is due to the non-zero contribution from F(a). 
In fact, since F(a) depends quadratically on the momentum Ef, the reduced phase space of 
the 2+1 Palatini theory with non-zero cosmological constant A is not naturally a cotangent 
bundle over a reduced configuration space. Thus, the result of Section 4 that the reduced 
phase space of the Chern-Simons theory based on G equals the reduced configuration space 
of the 2+1 Palatini theory based on the same G does not extend in general to the case A ^ 0. 

Nevertheless, we can still evaluate the Poisson brackets of the constraint functions F(a) 
and G(v). As in the A = case, we find that 

[G(y),G(w)} = G([v,w]), (5-18) 

where [v, w] 1 = C 1 jkv j w k is the Lie bracket of v 1 and w 1 , so v 1 i— ► G(v) is a representation 
of the Lie algebra Lq. Although F(a) has changed, we again find that 

{G(v),F(a)} = -F({v,a}), (5.19) 

where {v, a}i = C K jjv J aK is the co-adjoint bracket of v 1 and olj. However, the Poisson 
bracket of F(a) with F(f3) is no longer zero; it equals 

{F(a),F(P)} = -AG(e(a,P)), (5.20) 

where e(a,/3) 7 := e IJK ajPx- Thus, the totality of constraint functions is closed under 
Poisson bracket — i.e., they form a 1st class set. In fact, since (5.18), (5.19), and (5.20) do 
not involve any structure functions, the constraint functions form a Lie algebra with respect 
to Poisson bracket. The mapping (a,v) £ C* G x C G \— > (F(a),G(v)) is a representation of 
the Lie algebra £ag of the A- deformation, AG, of the Lie group G!^\ The F(a)'s play the 
role of "boosts" while the G(v)'s play the role of "rotations" in the A-deformation of G. 

5.2 Relationship to Chern-Simons theory 

In a manner similar to that used in subsection 4.4, we will now show that if G is any Lie 
group which admits an invariant, totally anti-symmetric tensor e IJK , then the 2+1 Palatini 
theory based on G with cosmological constant A is equivalent to Chern-Simons theory based 
on the A-deformation, AG, of the Lie group G. The actions for these two theories are the 
same modulo a surface term that does not affect the equations of motion. 

20 We will discuss the construction of the A-deformation of a Lie group G in the following section where 
we show the equivalence between the 2+1 Palatini theory based on G with cosmological constant A and 
Chern-Simons theory based on AG. When A = 0, AG is just the inhomogeneous Lie group IG constructed 
in subsection 4.4. 
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Given a Lie group G with an invariant, totally anti-symmetric tensor e , we first 
construct the A- deformation, AG, of G as follows: Form the direct sum vector space C G @Ca 
(having typical elements (a, v) 1 := (aj, v 1 ) and (f3, w) 1 := ((3j, w 1 )) and then define a bracket 
on C* G © C G via 

[(a, v), (P, w)]* := (-{v, P} + {w, a}, [v, w] - Ae(a, (3))\ (5.21) 

where [v,wY := C I JK v J w K , {v,(3}i := C K jiV J p K , and e(a, := e IJK aj(3 K - By inspec- 
tion, (5.21) is linear and anti-symmetric. By using the Jacobi identity C M \ijC n ' k]m = 
on Cq together with the anti-symmetry and invariance of e IJK , one can show that (5.21) 
satisfies the Jacobi identity as well. Thus, the vector space £ag := £-h © together with 
(5.21) is actually a Lie algebra. We call C\g the A- deformed Lie algebra associated with G. 
The A- deformation, AG, of G is obtained by exponentiating C\g- We can think of AG as 
an extension of the inhomogeneous Lie group IG in the sense that AG reduces to IG when 
A = 0. Note also that if G = SO{2, 1), then the above construction for AG yields 5*0(3, 1) 
if A < and 50(2,2) if A > 0. 

In addition to the above Lie algebra structure, C g @Cg is also equipped with a (natural) 
invariant, non-degenerate bilinear form 

kijia^Yip^y := arw 1 + (3^. (5.22) 

This is the same kij that we had when A = 0. As before, the existence of kij will allow us to 
construct Chern-Simons theory for AG. 

Given these remarks, we are now ready to verify that the 2+1 Palatini theory based on 
G with cosmological constant A is equivalent to Chern-Simons theory based on AG. Recall 
that 

%( 3 e, 3 A) = \ [ ^ 3 e a i( 3 FL - ±e" K \j 3 e cK ), (5.23) 

where 3 .A^ and \ a i are Cq- and /^-valued 1-forms. If we now construct the A-deformed Lie 
algebra Cag associated with G and define an /^AG-valued connection 1-form 3 A l a via 

:= ( 3 e a/ , 3 A{), (5.24) 

then a straightforward calculation along the lines of that used in subsection 4.4 shows that 
the Chern-Simons action 

AG S CS ( 3 A) = 1 -J m f] abc kij ( ^d h 3 Ai + I WfA, 34 C P) (5.25) 

equals 5a ( 3 e, 3 A) modulo a surface term which does not affect the Euler-Lagrange equations 
of motion. Specializing to the case G = 50(2, 1), we see that 2+1 gravity coupled to the 
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cosmological constant A is equivalent to Chern-Simons theory based on SO(3, 1) if A < or 
SO(2, 2) if A > 0. This was the observation of Witten @. 

5.3 2+1 Palatini theory coupled to a massless scalar field 

So far, we have seen that the 2+1 Palatini theory (with or without a cosmological constant 
A) is well-defined for a wide class of Lie groups. If A = 0, the Lie group G can be completely 
arbitrary; if A ^ 0, then G has to admit an invariant, totally anti-symmetric tensor e IJK . We 
are not forced to restrict ourselves to G = SO(2, 1). However, in order to couple fundamental 
matter fields with local degrees of freedom to 2+1 gravity via the 2+1 Palatini action, we 
will need to take G = 5*0(2,1). The matter actions require the existence of a spacetime 
metric g ab , and, as such, V a must have the interpretation of a co-triad. We will only consider 
coupling a massless scalar field to 2+1 gravity in this section — a similar treatment would 
work for Yang-Mills and Dirac fields as well. 

Let us first recall that the theory of a massless scalar field can be defined in n+1 
dimensions. If g ab denotes the inverse of the spacetime metric g ab , then the Klein-Gordon 
action SKG(g ab , 0) is defined by 

S KG (g ab A) ■= -Stt / ^=0 g ah d a <j>d h <l>, (5.26) 

where d a (j) denotes the gradient of 0. To couple the scalar field to gravity, we simply add the 
Klein-Gordon action (5.26) to the standard Einstein-Hilbert action 

S EH {g ab ) = I V^gR- (5.27) 

JM 

The total action ST{g ab , 0) is then given by the sum 

S T (g a \ 0) := S EH (g ab ) + S KG {g a \ 0), (5.28) 

and the Euler-Lagrange equations of motion are obtained by varying ST{g ab , 0) with respect 
to both g ab and 0. The variation of yields 

2 afe V a V fe = 0, (5.29) 

while the variation of g ab yields 

G ab = 87tT ab (KG). (5.30) 

V a is the unique, torsion-free spacetime derivative operator compatible with the metric g a b, 
and T a b(KG) is the stress-energy tensor of the massless scalar field. In terms of g ab and 0, 
we have 

T ab (KG) := d a <pd h <P - ^g ab d c <pd c <p. (5.31) 
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Now let us restrict ourselves to 2+1 dimensions and rewrite the above actions and equa- 
tions of motion in 2+1 Palatini form. As we saw in Section 3, the 2+1 Palatini action can 
be written as 

S P (% 3 A):=U f, ahc 3 e aI 3 F/ C , (5.32) 

I JM 

where 3 e a i is the co-triad related to the spacetime metric g a b via g a b :— 3 e I3 e J h r}u and 3 F^ b = 
2d[ a 3 Aq + e 1 JK 3 A^ 3 A^ is the internal curvature tensor of the generalized derivative operator 
^Da defined by 3 A I a . The Klein-Gordon action, viewed as a functional of 3 e^ and the scalar 
field 0, is given by 

SkgCc^) = -8tt / g ab d a (t>d h <t>. (5.33) 

JM 

Note that although the Klein-Gordon action depends on the co-triad 3 e aI through its depen- 
dence on y/^^g and g ab , it is independent of the connection 1-form 3 A^. In fact, of all the 
fundamental matter couplings, only the action for the Dirac field would depend on 3 A T a . 
Given (5.32) and (5.33), we define the total action as the sum 

S T (% 3 A,<j>) ■= S P ( 3 e, 3 A) + \s KG {\<P). (5.34) 

The additional factor of 1/2 is needed in front of Skg^i 0) so that the above definition 
of the total action will be consistent with the definition of Sp( 3 e, 3 A). The Euler-Lagrange 
equations of motion are obtained by varying Sr( 3 e, 3 A, 0) with respect to each field. Varying 
gives 

S afe V a V b = 0, (5.35) 

while varying 3 A ! a and 3 e a j imply 

*D b (v abc 3 e c /) = and (5.36) 

~abc 3^/ _ 8n ^j(3 e al g bc _ 2 3 £ bl g acj Q^Q^ = Q ? (5.37) 

respectively. Note that equation (5.35) is just the standard equation of motion for 0, while 
equation (5.36) implies that 3 A^ = T^, as in the vacuum case. Substituting this result for 
3 A T a back into (5.37) and contracting with 3 ef gives 

Q ad = 87rT ad( KG y ( 5>38 ) 

These are the desired results. 

To put this theory in Hamiltonian form, we will basically proceed as we have in the past, 
but use additional structure provided by the spacetime metric g a b- We will assume that 
M is topologically E x R for some spacelike submanifold E and assume that there exists 
a time function t (with nowhere vanishing gradient (dt) a ) such that each t = const surface 
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S t is diffeomorphic to E. t a will denote the time flow vector field (t a (dt) a = 1), while n a 
will denote the unit covariant normal to the t = const surfaces. n a := g ab n b will be the 
associated future-pointing timelike vector field (n a n a = —1). Given n a and n a , it follows 
that q b := 5 b + n a n b is a projection operator into the t = const surfaces. We can then define 
the induced metric q ab , the lapse N, and shift N a as we did for the standard Einstein-Hilbert 
theory in Section 2. Recall that t a = Nn a + N a , with N a n a = 0. 

Now let us write the total action (5.34) in 2+1 form by decomposing each of its pieces. 
Using rj abc = Zforf^dt, it follows that 

S P ( 3 e, 3 A) = JdtJ^ i( 3 e • t)/?f b Fi, + E T C r A l a - E,V a { 3 A ■ t) 1 , (5.39) 

where ( 3 e •*)/:= t a 3 e aI , Ef := rj ab 3 e bI , ( 3 A ■ t) 1 := t a 3 A J a , and A{ := q b a 3 A{. To obtain (39), 
we used the fact that £ t - q b = 0. Note also that F^ b := q c a q b 3 F/ d is the curvature tensor of 
the generalized derivative operator V a (:— q b ^Dj,) on E associated with A T a . Since 

l -{ 3 e ■ t)jv ab F^ b = —^Ne IJK EjEjFabx - N a E b F^ b (5.40) 

(where N := q~^N), we see that (modulo a surface integral) the Lagrangian Lp of the 2+1 
Palatini theory is given by 

Lp = L -\^ IJK Z a iE b jF abK - N a E b F' b + Ef£ t A{ + (V a Ef) ( 3 A ■ t) 1 . (5.41) 

Similarly, using g ab = q ab — n a n b and the decomposition \f—g = Ny/qdt, it follows that 

Skg{\4>) = -8tt J dt ^{Nr b d a( pd b 0-N-\£ r 0-N a d a( p) 2 }, (5.42) 

where q ab := qq ab (= E I E bI ). Thus, the Klein-Gordon Lagrangian L KG is simply given by 

Lkg = -Stt / {Nr b d a <f>d b <f> - %-\C;4>- N a d a <P) 2 }. (5.43) 

The total Lagrangian L T is the sum L T = L P + \L KG and is to be viewed as a functional of 
the configuration variables ( 3 A ■ t) 1 , AJ", JV°, and their first time derivatives. 

Following the standard Dirac constraint analysis, we find that 

tt := = STriV-^A- - iV^ a 0) (5.44) 



is the momentum canonically conjugate to 0. Since this equation can be inverted to give 

1 

8tt 



= —Nn + N a d a 4>, (5.45) 
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it does not define a constraint. On the other hand, Ef is constrained to be the momentum 
canonically conjugate to A^, while ( 3 A -t) 1 , N, and N a play the role of Lagrange multipliers. 
The resulting total phase space (T T ,f2 T ) is coordinatized by the pairs of fields (A T a ,Ej) and 
(0, 7f ) with symplectic structure 

tt T = J &E? A &A{ + Tr(dLr A &<f>). (5.46) 

The Hamiltonian is given by 

H T (A, E, 0, E) = ( N (^-e IJK EjEjF a bK + (4^.0^0 + -I_ n 2 )) 

js y z 107T ' (5.47) 

+ N a (E b X b + ^0) - (V a Ef)( 3 A ■ t) 1 , 
We shall see that this is just a sum of 1st class constraint functions associated with 

±e IJK E*E b jF abK + (4vRfU<W + -L n 2 ) » 0, (5.48) 

Z l07T 

+ vf9 a w 0, and (5.49) 
V a E a j w 0. (5.50) 

These are the constraint equations associated with the Lagrange multipliers N, N a , ( 3 A-tY, 
respectively. 

Two remarks are in order: First, note that just as we found for the 2+1 Palatini theory 
with or without a cosmological constant A, the constraint equations for the 2+1 Palatini 
theory coupled to a Klein-Gordon field are polynomial in the basic canonically conjugate 
variables (A^Ef) and (0, n). Since the Hamiltonian is just a sum of these constraints, it 
follows that the evolution equations will be polynomial as well. Second, since the constraint 
equations do not involve the inverse of Ej, the above Hamiltonian formulation is well-defined 
even if Ef is non-invertible. Thus, we have a slight extension of the standard 2+1 theory of 
gravity coupled to a massless scalar field. It can handle those cases where the spatial metric 
~ab ._ fiafibi b ecome s degenerate. 

Our next goal is to verify the claim that the constraint functions associated with (5.48)- 
(5.50) form a 1st class set. To do this, we let v 1 (which takes values in the Lie algebra of 
5*0(2, 1)), N, and iV a be arbitrary test fields on E. Then we define 

C(N) := N(^e IJK EjEjF a bK + (4^40^0+^ ^f 2 )), (5.51) 
C\N) := f N^E'F^ + W^), and (5.52) 
G(v) := / v\V a E% (5.53) 



s 
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to be the scalar, vector, and Gauss constraint functions. 

As we saw in subsection 3.3 for the 2+1 Palatini theory, the subset of Gauss constraint 
functions form a Lie algebra with respect to Poisson bracket. Given G(v) and G(w), we have 

{G(v),G(w)} = G([v,w]), (5.54) 

where [v, w} 1 is the Lie bracket in £50(2,1)- Thus, the mapping v 1— »■ G(v) is a representation 
of the Lie algebra £50(2,1)- Given its geometrical interpretation as the generator of internal 
rotations, it follows that 

{G(v),C(N)} = and (5.55) 
{G(v),C'(N)} = 0, (5.56) 

as well. 

Since one can show that the vector constraint function does not by itself have any direct 
geometrical interpretation (see, e.g., [0), we will define a new constraint function C(N) by 
taking a linear combination of the vector and Gauss constraints. We define 

C(N) ■= C'{N) - G(N), (5.57) 

where N 1 := N a A I a . We will call C(N) the diffeomorphism constraint function since the 
motion it generates on phase space corresponds to the 1-parameter family of diffeomorphisms 
on £ associated with the vector field N a . To see this, we can write 

C(N) : = C\N) - G{N) 

N'i&jF^ + ndat)- [ N^V^) 

(5 58) 

N a (E b jFi b - AlV b E b T ) + m a d a <P 



where the Lie derivative with respect to iV a treats fields having only internal indices as 
scalars. To obtain the last line of (5.58), we ignored a surface integral (which would vanish 
anyways for iV a satisfying the appropriate boundary conditions). By inspection, it follows 
that A{ 1— >■ + eCfiAa + 0(e 2 ), etc. Using this geometric interpretation of C(N), it follows 
that 

{C(N),G(v)} = G(C$v), (5.59) 
{C(N),C(M)} = C(C^M), and (5.60) 
{C(N),C(M)} = C([N,M]). (5.61) 
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We are left to evaluate the Poisson bracket {C(N),C(M)} of two scalar constraints. 
After a fairly long but straightforward calculation, one can show that 



{C(N), C(M)} = C\K) (= C{K) + G(K, K)) , 



(5.62) 



where K a := q ab (Nd b M - Md b N) and q ab = EfE bI . This result makes crucial use of the 



This is a property of the structure constants e 1 jk '■= £ IMN VmjVnk of the Lie algebra of 
5*0(2, 1). Thus, the constraint functions are closed under Poisson bracket — i.e., they form 
a 1st class set. Note, however, that since the vector field K a depends on the phase space 
variable Ej, the Poisson bracket (5.62) involves structure functions. The constraint functions 
do not form a Lie algebra. This result is similar to what we found for the standard Einstein- 
Hilbert theory in subsection 2.3. 

It is interesting to note that even if we did not couple matter to the 2+1 Palatini theory, 
but performed the Legendre transform as we did above (i.e., using the additional structure 
provided by the spacetime metric g ab : = V a 3 e^/j), we would still obtain the same Poisson 
bracket algebra. The constraint functions would still fail to form a Lie algebra due to the 
structure functions in (5.62). At first, something seems to be wrong with this statement, 
since we saw in subsection 3.3 that the constraint functions G(v) and F(a) of the 2+1 
Palatini theory form a Lie algebra with respect to Poisson bracket. One may ask why the 
constraints functions obtained via one decomposition of the 2+1 Palatini theory form a Lie 
algebra, while those obtained from another decomposition do not. 

The answer to this question is actually fairly simple. Namely, it is easy to destroy the 
"Lie algebra- ness" of a set of constraint functions. If & (i — 1, • • • , m) denote m constraint 
functions which form a Lie algebra under Poisson bracket (i.e., {(pi, <pj} = C-ijipk where C-jj 
are constants), then a linear combination of these constraints, \%_ = ^jr-'Pji wm n °t m general 
form a Lie algebra if A^- are not constants on the phase space. In essence, this is what 
happens when one passes from the G(v ) and F(a) constraint functions of subsection 3.3 to 
the G(v ), C(N), and C(N) constraint functions of this subsection. The transition from F(a) 
to C(N) and C(N) involve functions of the phase space variables. 

6. 3+1 Palatini theory 

In this section, we will describe the 3+1 Palatini theory. In subsection 6.1, we define the 
3+1 Palatini action and show that the Euler-Lagrange equations of motion are equivalent 



fact that 






(5.63) 
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to the standard vacuum Einstein's equation. In subsection 6.2, we will follow the standard 
Dirac constraint analysis to put the 3+1 Palatini theory in Hamiltonian form. We obtain a 
set of constraint equations which include a 2nd class pair. By solving this pair, we find that 
the remaining (1st class) constraints become non-polynomial in the (reduced) phase space 
variables. In essence, we are forced into using the standard geometrodynamical variables of 
general relativity. In fact, as we shall see in subsection 6.3, the Hamiltonian formulation of 
the 3+1 Palatini theory is just that of the standard Einstein-Hilbert theory. Thus, the 3+1 
Palatini theory does not give us a connection-dynamic description of 3+1 gravity. 

Much of the material in subsection 6.2 is based on an an analysis of the 3+1 Palatini 
theory given in Chapter 4 of [Q. 

6.1 Euler-Lagrange equations of motion 

To obtain the Palatini action for 3+1 gravity, we will first write the standard Einstein- 
Hilbert action 

S EH (g ab ) = I V^R (6.1) 

in tetrad notation. Using 

RabI J = Rabc 1 Gj Is^ (6.2) 

(which relates the internal and spacetime curvature tensors of the unique, torsion-free gen- 
eralized derivative operator V a compatible with the tetrad %") and 

tabcd — e a e b e c e d e IJKL yo.6) 

(which relates the volume element e a b c d of g a b = ^„ ^bVu to the volume element eijKL of 
77/ j), we find that 

V=gR = \v abcd eiJKL *e{ V R cd KL . (6.4) 

Thus, viewed as a functional of the co-tetrad the standard Einstein-Hilbert action is 
given by 

SehCc) = \ £ fj abcd eiJKL H A e J b R cd KL . (6.5) 

To obtain the 3+1 Palatini action, we simply replace R a bi J in (5) with the internal 
curvature tensor 4 F a bj J of an arbitrary generalized derivative operator defined by 

*D a k! := d a kj + 4 A aI J kj. (6.6) 

We define the 3+1 Palatini action to be 

S P (% M) := I / A;« K H 4 F cd KL , (6-7) 
8 Jm 
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where 4 F abI J = 2d[ a 4 A b ] I J + [ 4 A a , 4 A b ]j J . Just as we did for the 2+1 Palatini theory, we have 
included an additional factor of 1/2 in definition (6.7) so our canonically conjugate variables 
will agree with those used in the literature (see, e.g., Note also, that as defined above, 
knows how to act only on internal indices. We do not require that know how to 
act on spacetime indices. However, we will often find it convenient to consider a torsion-free 
extension of ^Va to spacetime tensor fields. All calculations and all results will be independent 
of this choice of extension. 

Since the 3+1 Palatini action is a functional of both a co-tetrad and a connection 1-form, 
there will be two Euler-Lagrange equations of motion. When we vary Sp( 4 e, 4 A) with respect 
to 4 e^ and 4 A a IJ , we find 

fVi^V^O and (6.8) 

4 ef = 0, (6.9) 

respectively. The last equation requires a torsion-free extension of ^Da to spacetime tensor 
fields, but since the left hand side of (6.9) is the divergence of a skew spacetime tensor 
density of weight +1 on M, it is independent of this choice. Noting that rj abcd euKL 4 ef 4 e L d = 
4( 4 e) 4 e^ 4 e b j (where ( 4 e) := y/—g), we can rewrite (6.9) as 

%> (( 4 e) 4 e? 4 e?) = 0. (6.10) 

This equation is identical in form to equation (3.18) obtained in Section 3 for the 2+1 Palatini 
theory. 

By following exactly the same argument used in subsection 3.1, equation (6.10) implies 
that 4 A aI J = T aI J , where T aI J is the internal Christoffel symbol of V a . Using this result, the 
remaining Euler-Lagrange equation of motion (6.8) becomes 

V abcd e IJK L 4 ei Rcd KL = 0. (6.11) 

When (6.11) is contracted with 4 e eI , we get G ae = 0. Thus, we can produce the 3+1 vacuum 
Einstein's equation starting from the 3+1 Palatini action given by (6.7). Note that just as 
in the 2+1 theory, the equation of motion (6.9) for A a j J can be solved uniquely for 4 A a j J 
in terms of the remaining basic variables The pulled-back action S_p( 4 e) defined on the 
solution space 4 A a i J = T a i J is just 1/2 times the standard Einstein-Hilbert action Seh(S) 
given by (6.5). 

6.2 Legendre transform 

To put the 3+1 Palatini theory in Hamiltonian form, we will use the additional structure 
provided by the spacetime metric g ab . We will assume that M is topologically E x R for 
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some spacelike submanifold £ and assume that there exists a time function t (with nowhere 
vanishing gradient (dt) a ) such that each t = const surface S t is diffeomorphic to E. t a 
will denote the time flow vector field (t a (dt) a = 1), while n a will denote the unit covariant 
normal to the t = const surfaces. n a := g ab nb will be the associated future-pointing timelike 
vector field (n a n a = —1). Given n a and n a , it follows that q b := 5 b + n a n b is a projection 
operator into the t = const surfaces. We can then define the induced metric q a b, the lapse 
N, and shift N a as we did for the standard Einstein-Hilbert theory in Section 2. Recall that 

t a = Nn a + N a^ with jya^ = Q 

Now let us write (6.7) in 3+1 form. Using the decomposition fj abcd = At^rf'^dt (where 
~abc j g L ev i-Civita tensor density of weight +1 on £), we get 



C 4j\ _ 1 /" -abcd^ AJ A J Arp KL 

bp{ e, A) — - i] e IJKL e a e b b cd 

o JM 

= l dt L ~S e ' t )' € "^ e hFcd KL + \&ijC* K 1J - \E a i.jV a ^A ■ t) IJ , 

(6.12) 

where (% • i)j := T % aJ , £ a 7J := | e/Ji rf fec %f ^, ( 4 A ■ t) IJ := t a 4 A a IJ , A a IJ := q b 4 A b IJ , 
and e T a := q b a 4 e T b . To obtain the last line of (6.12), we used the fact that £ t - q b = 0. Note 
also that F ab IJ := q c a q b A F c d IJ is the curvature tensor of the generalized derivative operator 
('■— Qa *Db) on E associated with A a IJ . Since 

\Ce ■ tYe IJKL fj bcd e J b F cd KL = -±NTr(E a E b F ab ) + ^N a Tr(E b F ab ) (6.13) 

(where N := q~^N and Tr denotes the trace operation on internal indices), we see that 
(modulo a surface integral) the Lagrangian L P of the 3+1 Palatini theory is given by 

L p = [ -\NTi{E a E b F ab ) + lN a Tr(E b F ab ) 

* ( 614 ) 

+ -E a u C r A a IJ + -(P a ^„)(t4 • t) IJ . 

The configuration variables of the theory are ( A A ■ t) IJ , N, N a , A a IJ , and E a jj. 

But before we perform the Legendre transform, we should note that the configuration 
variable E a jj is not free to take on arbitrary values. In fact, from its definition 

E a u := \euKLV abc 4 ef 4 e c L , (6.15) 

one can show that 

4> ab := e IJKL E a uE b KL = and Tr(£ a £ b ) > 0. (6.16) 
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The second condition follows from the fact that Tr(E a E b ) = 2q ab (= 2qq ab ), where q ah is 
the inverse of the induced positive-definite metric q a b on E. Thus, the starting point for the 
Legendre transform is L P togetherwith the primary constraint <p ab = 0. Since the inequality is 
a non-holonomic constraint, it will not reduce the number of phase space degrees of freedom. 

If we now follow the standard Dirac constraint analysis, we find that ( 4 A ■ t) IJ , N, and 
N a play the role of Lagrange multipliers. Their associated constraint equations (which arise 
as secondary constraints in the analysis) are 

Tr(E a E b F ab ) « 0, (6.17) 
Tr(E b F ab ) « 0, and (6.18) 
V a E a u « 0. (6.19) 

There is also a primary constraint which says that \E a u is the momentum canonically 
conjugate to A a IJ . By demanding that the Poisson bracket of this constraint with the total 
Hamiltonian and the Poisson bracket of (f> ab with the total Hamiltonian be weakly zero, we 
find that 

X ab ■= e IJKL {V c E a u ) [E b , E c ]kl + (a^b)^0. (6.20) 

This is an additional secondary constraint which must be appended to constraint equations 
(6.16)-(6.19). In virtue of <p ab = 0, the expression for x ab is independent of the choice 
of torsion-free extension of V a to spacetime tensor fields. If we further demand that the 
Poisson bracket of x ab with the total Hamiltonian be weakly zero, we find nothing new — i.e., 
there are no tertiary constraints. 

Let us summarize the situation so far: Out of the original set of configuration variables 
{( 4 A ■ t) IJ , N, N a , A a IJ , E a u}, the first three are non-dynamical. We also found that 
\E a u is the momentum canonically conjugate to A a IJ . Thus, the phase space (T' P , il' P ) of 
the 3+1 Palatini theory is coordinatized by the pair (A a IJ ,E a jj) and has the symplectic 
structure 

n' P = I [ &E a u A dA a IJ . (6.21) 

The Hamiltonian is given by 



H P (A, E) = jf l -NTv{E a E b F ab ) - l -N a Tv{E b F ab ) - \{V a E a U )( 4 A ■ t) 1 



(6.22) 



+ ^€ rjKL E a IJ E b 



IJ-n KLi 



where Apt is a Lagrange multiplier. The constraints of the theory are (6.16)-(6.20). Note also 
that at this stage of the Dirac constraint analysis all constraint (and evolution) equations 
are polynomial in the canonically conjugate variables. 
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The next step in the Dirac constraint analysis is to evaluate the Poisson brackets of the 
constraints and solve all 2nd class pairs. Since only 

& b (x), X cd (y)}^0, (6.23) 

we just have to solve constraint equations (6.16) and (6.20). As shown in Chapter 4 of [^J, 
the most general solution to (6.16) is 

E a 7J =:2E [ > J] , (6.24) 

for some unit timelike covariant normal nj {njn 1 = — 1) with Ef invertible and Efn 1 = 0. 
By writing E a jj in this form, we see that the original 18 degrees of freedom per space point 
for E a u has been reduced to 12. Note also that EfE bI = q ab , so Ef is in fact a (densitized) 
triad. 

Given (6.24), the most convenient way to solve (6.20) is to gauge fix the internal vector 
n 1 . This will further reduce the number of degrees of freedom of E a jj to 9, since now 
only Ef will be arbitrary. However, gauge fixing n 1 requires us to solve the boost part of 

(6.19) relative to ni as well.f^] We can only keep that part of (6.19) which generates internal 
rotations leaving n 1 invariant. 

To solve these constraints, let us first define an internal connection 1-form K a j J via 

V a kj =: D a kj + K aI J kj, (6.25) 

where D a is the unique, torsion-free generalized derivative operator on £ compatible with 
the (densitized) triad Ef and the gauge fixed internal vector n 1 . Then constraint equation 

(6.20) and the boost part of (6.19) become 

x ab = -4e IJK K cI L E ( L a E b ]E c K « and (6.26) 
(V a E a u)n J = -K aM N E a N qf « 0, (6.27) 

where qj := 5j + n 7 nj. By using the invertibility of the (densitized) triad Ef, one can then 
show (again, see Chapter 4 of 0) that (6.26) and (6.27) imply that K a IJ also be pure boost 
with respect to ni — i.e., that K a IJ have the form 

K a " -. 2K^n J \ (6.28) 

with K^nj = 0. Since D a is determined completely by Ef and rij, the original 18 degrees of 
freedom for A a IJ has also been reduced to 9 degrees of freedom per space point. The infor- 
mation contained in A a IJ (which is independent of Ef and rii) is completely characterized 

21 The boost part of any anti-symmetric tensor Ajj relative to ri] is defined to be Ajjn J . 
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by K*. To emphasize the fact that Efn 1 = and K^ni = 0, we will use a 3- dimensional 
abstract internal index i and write Ef and K l a in what follows. 

Thus, after eliminating the 2nd class constraints, the phase space (Tp,flp) of the 3+1 
Palatini theory is coordinatized by the pair (Ef, K^) and has the symplectic structure 



n P = / &KI A &Ef. (6.29) 

The Hamiltonian is given by 

H P (E, K) = J I#( - q K-2E$\K\Kt) - 2N a E\D {a K\ 



(6.30) 



+ ( 4 A ■ tfEZK, 



where 71 denotes the scalar curvature of D a . This is just a sum of the 1st class constraints 
functions associated with 



C(E, K) := -qTZ - 2E?MK\Kl « 0, (6.31) 



3 



C a (E, K) := ±E\D {a K\ w 0, and (6.32) 
G ij (E,K):=-Ef i K aj] ^0. (6.33) 

(The overall numerical factors have been chosen in order to facilitate the comparison with 
the standard Einstein-Hilbert theory.) Note that constraint equations (6.31), (6.32), and 
(6.33) are the remaining constraints (6.17), (6.18), and the rotation part of (6.19) relative 
to rii expressed in terms of the phase space variables (Ef,K l a )!^ We will call (6.31), (6.32), 
and (6.33) the scalar, vector, and Gauss constraints for the 3+1 Palatini theory. 

Note that as a consequence of eliminating the 2nd class constraints by solving (6.16), 
(6.20), and the boost part of (6.19) relative to nj, the constraint equations (6.31)-(6.33) 
(and hence the evolution equations generated by the Hamiltonian) are now non-polynomial 
in the canonically conjugate pair (Ef, K^). This is due to the dependence of 1Z on the inverse 
of Ef. In fact, since Ef must be invertible, we are forced to take Ef as the configuration 
variable of the theory. We are led back to a geometrodynamical description of 3+1 gravity. 
Thus, the Hamiltonian formulation of the 3+1 Palatini theory has the same drawback as the 
Hamiltonian formulation of standard Einstein-Hilbert theory. As we shall see in the next 
subsection, these theories are effectively the same. 

6.3 Relationship to the Einstein-Hilbert theory 



22 The rotation part of any anti-symmetric tensor Au relative to ni is given by qf 1 Amn where qj :— 
5 j + "nJnj. 
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In this subsection, we will not explicitly evaluate the Poisson bracket algebra of the 
constraint functions for the 3+1 Palatini theory. Rather, we will describe the relationship 
between the constraint equations (6.31)-(6.33) and those of the standard Einstein-Hilbert 
theory. We shall see that if we solve the first class constraint (6.33) by passing to a reduced 
phase space, we recover the Hamiltonian formulation of the standard Einstein-Hilbert theory 
in terms of the induced metric q ab and its canonically conjugate momentum p ab . 

To do this, let us first define a tensor field M ab (of density weight -1 on S) via 

M ab := KlE bl) (6.34) 

where E\ is the inverse of the (densitized) triad Ef. Then in terms of M ab) one can show 
that constraint equation (6.33) is equivalent to 

M [ab] « 0. (6.35) 

Thus, the constraint surface in Tp defined by (6.33) will be coordinatized in part by the 
symmetric part of M ab — i.e., by K ab := M {ab) . 

But we are not yet finished. Since (6.33) is a 1st class constraint, we must also factor- 
out the constraint surface by the orbits of the Hamiltonian vector field associated with the 
constraint function 

'■— Gij(E, K)A ij (= j^-E%K aj Wy (6.36) 

(Here A JJI = A^l denotes an arbitrary anti-symmetric test field on E.) Since it is fairly 
easy to show that G(A) generates (gauge) rotations of the internal indices (i.e., E°; h- > + 
eAi j Ej+0(e 2 ) and K l a i-> K l a — eA/i^ + 0(e 2 )), the factor space will be coordinatized by K ab 
and the gauge invariant information contained in E?. This is precisely q ab = E^E bt . Thus, 
the reduced phase space (fp,Qp) is coordinatized by the pair (q ab ,K ab ) and has symplectic 
structure 



ttp = -J^K ab A dig . (6.37) 

All we must do now is make contact with the usual canonical variables of the standard 
Einstein-Hilbert theory. To do this, let us work with the undensitized fields q ab and K ab: 
and lower and raise their indices, respectively. Then in terms of p ab defined by 

p ab ■= ^(K ab - Kq ab ), (6.38) 

one can show that 



= f dp ab Adq ab (6.39) 
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and 



C(q, p) = -qK + (p ab p ab ~ ^f) « 0, (6.40) 
C a (q,p) = -2q ab D c p bc ^V. (6.41) 

Up to overall factors, these are just the symplectic structure ^Ieh and scalar and vector 
constraint equations of the standard Einstein-Hilbert theory described in Section 2. The 
factor of —1/2 in the symplectic structure is due to the combination of using an action 
which is 1/2 the standard Einstein-Hilbert action and using a (densitized) triad Ef instead 
of a covariant metric q ab as our basic dynamical variables. Thus, we see that the Hamiltonian 
formulation of the 3+1 Palatini theory is nothing more than the familiar geometrodynamical 
description of general relativity. 



7. Self-dual theory 

In this section, we will describe the self-dual theory of 3+1 gravity. This theory is similar 
in form to the 3+1 Palatini theory described in the previous section; however, it uses a 
self- dual connection 1-form as one of its basic variables. We define the self-dual action for 
complex 3+1 gravity and show that we still recover the standard vacuum Einstein's equation 
even though we are using only half of a Lorentz connection. We then perform a Legendre 
transform to put the theory in Hamiltonian form. In terms of the resulting complex phase 
space variables, all equations of the theory are polynomial. This simplification gives the 
self-dual theory a major advantage over the 3+1 Palatini theory. As noted in the previous 
section, the Hamiltonian formulation of the 3+1 Palatini theory reduces to that of the 
standard Einstein-Hilbert theory with its troublesome non-polynomial constraints. 

As mentioned in footnote 3 in Section 1, to obtain the phase space variables for the real 
theory, we must impose reality conditions to select a real section of the original complex 
phase space. At the end of subsection 7.2 we will describe these conditions and discuss 
how they are implemented. The need to use reality conditions is a necessary consequence 
of using an action principle to obtain the new variables for real 3+1 gravity. Although we 
mention here that it is possible to stay within the confines of the real theory by performing 
a canonical transformation on the standard phase space of real general relativity, we will not 
follow that approach in this paper. (Interested readers should see JTJ for a detailed discussion 
of Ashtekar's original approach.) Rather, we will start with an action for the complex theory 
and obtain the new variables as outlined above. Henceforth, the co-tetrad 4 e^ will be assumed 
to be complex unless explicitly stated otherwise. 
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7.1 Euler-Lagrange equations of motion 

To write down the self-dual action for complex 3+1 gravity, all we have to do is replace 
the (Lorentz) connection 1-form 4 A a j J of the 3+1 Palatini theory by a self-dual connection 
1-form ^ 4 A aI J and let the co-tetrad *e a become complex. We define the self-dual action to be 

S SD (% «A) := \ I v abcd e IJK L V a % ^F cd KL , (7.1) 

where ~ H F abI J = 2d[ a +4 A b ] I J + \^A a , ^4 b ]/ J is the internal curvature tensor of the self-dual 
generalized derivative operator ~ i4 V a defined by 

*V a k! := d a kj + ^Ajkj. (7.2) 

Some remarks are in order: 

1. We will always take our spacetime manifold M to be a real 4-dimensional manifold. 
Complex tensors at a point p G M will take values in the appropriate tensor product 
of the complexified tangent and cotangent spaces to M at p. The fixed internal space 
will also be complexified; however, the internal Minkowski metric rju will remain real. 
Since the co-tetrad ^e l a are allowed to be complex, the spacetime metric g a b defined by 
9ab '■= 4 Cq 4 e b rju will also be complex. 

2. Although we can no longer talk about the signature of a complex metric g ab , com- 
patibility with a complex co-tetrad 4 e^ still defines a unique, torsion-free generalized 
derivative operator V„. Thus, the complex Einstein tensor G ab :— R a b — \Rg a b is 
well-defined; whence the complex equation of motion G a b = makes sense. It is this 
equation that defines for us the complex theory of 3+1 gravity. 

3. When we say that the connection 1-form +A A a i J (or any other generalized tensor field) 
is self-dual, we will always mean with respect to its internal indices. Thus, the notion 
of self-duality makes sense only in 3+1 dimensions and applies only to generalized 
tensor fields with a pair of skew-symmetric internal indices, T a "' b c ... dIJ = T a '" b c ... d [ IJ ]. 
The dual of T a '" b c ... dIJ , denoted by *T a '" b c ... dIJ , is defined to be 

c-dIJ ■— -£IJ 1 c-dKL, yi-o) 

where the internal indices of €ukl are raised with the internal metric n IJ . Since rju 
has signature ( — h ++), it follows that the square of the duality operator is minus 
the identity. Hence, our definition of self-duality involves the complex number %. We 
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define T a '" b c ... dIJ to be self- dual if and only ifQ 

Ta-b ■rria—b 
c-dIJ — " c-dIJ- 



(7.4) 



Thus, self-dual fields in Lorentzian 3+1 gravity are necessarily complex. 

4. Given any generalized tensor field T a '" b c ... dIJ = T a '" b c .- d[/j], we can always decompose 
it as 

Ta---b -+rpa---b . —rna---b r\ 

c-dIJ — 1 c-dIJ T 1 c-dIJ, l'-CJ 

where 

•V" 'V...,/,./ := \{T a - b c ... dIJ - iT a - b c ... dU ) and (7.6a) 

Since I 'V" 'V... (// 7 ) = i+T a - b c ... dIJ and ( •/•" '•,...,„./) = -i T a - b c ... dIJ , it follows that 
+ T a '" t 'c-du and ~T a "' b c-diJ are the self-dual and antt self-dual parts of T a "' b c ... d ij. 
Equations (7.6a) and (7.66) define the self-duality and anti self-duality operators + 
and _ . 

5. The generalized derivative operator ^D a is said to be self-dual only in the sense that it 
is defined in terms of a self-dual connection 1-form 44 A a / J . As in many of the previous 
theories, J ^D a (as defined by (7.2)) knows how to act only on internal indices. But 
as usual, we will often find it convenient to consider a torsion-free extension of JvA D a 
to spacetime tensor fields. All calculations and all results will be independent of this 
choice of extension. Note also that the internal curvature tensor of the generalized 
derivative operator J ^ r D a is given by 

*F M J = 2 + [*A a , ^A b )j J . (7.7) 

Since one can show that the (internal) commutator of two self-dual fields is also self- 
dual, it follows that ~^F a bi J is self-dual with respect to its internal indices. 

Given these general remarks, we are now ready to return to the self-dual action (7.1) 
and obtain the Euler-Lagrange equations of motion. Varying Ssd^, 44 ^4) with respect to 4 e ! a 
gives 

V^e IJKL % J b «F cd KL = 0, (7.8) 

23 T a - b c ... dIJ is defined to be anti self-dual if and only if *T a - b c ... d ij = 'T" "'\ ...,//./ . The choice of +i 
for self-dual and — i for anti self-dual is purely convention. I have chosen our conventions to agree with those 
of 1- 
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while varying Ssoi^^A) with respect to ^Aj J gives 

^(( 4 e) + ( 4 et a4 e3))=0. (7.9) 

To obtain (7.9), we used the fact that rj abcd euKL ^ ^ = 4( 4 e) 4 e^ 4 e) (where (1s) := y/—g). 
Note also that we are forced to take the self-dual part of A e^ 4 e] since the variations S^Aj 3 
are required to be self-dual. This is the distinguishing feature between the self-dual and 
3+1 Palatini equations of motion. Finally note that although (7.9) requires a torsion- free 
extension of ^T> a to spacetime tensor fields, it is independent of this choice since the left 
hand side is the divergence of a skew spacetime tensor density of weight +1 on M. 

We would now like to show that (7.9) implies that ~^V a is the self-dual part of the unique, 
torsion- free generalized derivative operator V a compatible with But since we are working 
with self-dual fields, the argument used for the 2+1 and 3+1 Palatini theories does not yet 
apply. We will have to do some preliminary work before we can use those results. 

If T aI J denotes the internal Christoffel symbol of V , we define the self-dual part of 
V a by 

* i Vaki:=d a k I + + r aI J kj, (7.10) 

where + T aI J is the self-dual part of T aI J . The difference between JrA T> a and ^\7 a is then 
characterized by a generalized tensor field +4 C a i J defined by 

^Dakj^.^kj + ^Ca/kj. (7.11) 

Note that +4 C a / J is self-dual as the notation suggests. In fact, +4 C a / J = ^A^ 3 — + r a / J . Now 
let us write (7.9) in terms of ~^V a and +4 C a / J . Using (7.11) to expand the left hand side of 
(7.9), we get 

M (( 4 e) + ( 4 e [ j VJ)) + ( 4 e) (^C bI K +( 4 e£ VJ) + «C„ K \^ 4 4)) = 0. (7.12) 

Since + V a ki = V a ki — ~T a i J kj (where ~T a j J is the anti self-dual part of the internal Christof- 
fel symbol r a / J ), and since the last two terms of (7.12) can be written as the (internal) 
commutator of ^Caij and + ( 4 ef 4 e b J), we get 

- -r b , +(V a yi)]^ + +(V a yi)]^ = o. (7.i3) 

The first commutator vanishes since ~T b jj is anti self-dual while + ( 4 e^' 4 e b }) is self-dual; in the 
second commutator, Vj!) can be replaced by %j Q VjL Thus, (7.13) reduces to 

^,(V a4 e fe ])] 7j = 0. (7.14) 
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This is exactly the form of the equation found in the 2+1 Palatini theory with +4 Cbi J replacing 
? Cbi J ■ (See equation (3.20).) We can now follow the argument given there to conclude that 
^C aI J = 0. Thus, ^A aI J = + T aI J as desired. 

By substituting the solution ^A aI J = + T aI J into the remaining equation of motion (7.8), 
we get 

V abcd euKL V b + R cd KL = (7.15) 



where + R c d is the self-dual part of the internal curvature tensor R c d of V a . Then by 
using the definition of + R c d KL , we see that equation (7.15) becomes 



r\—^zpbcd c 4J/p KL * KL d MN\ 

u — -rj euKL e b {K c d — -e MN-ticd ) 

1 1 (7.16) 

_ -L-abccT A J p KL 

— -jV tlJKL %ticd , 



where the second term on the first line vanishes by the Bianchi identity R\ a bc]d — 0. When 
(7.16) is contracted with 4 e e/ , we get G ae = 0. Thus, the self-dual action (7.1) reproduces 
the vacuum Einstein's equation for complex 3+1 gravity. 

Since this is an important — yet somewhat suprising — result, it is perhaps worthwhile to 
repeat the above argument from a slightly different perspective. First note that the self-dual 
action (7.1) and the 3+1 Palatini action (6.7) differ by a term involving the dual of the 
curvature tensor A F c d KL . This extra term in the self-dual action is not a total divergence and 
thus gives rise to an additional equation of motion that is not present in the 3+1 Palatini 
theory. This equation of motion also involves the dual of the curvature tensor. (Compare 
equations (7.8) and (6.8).) However, as we showed above, if we solve (7.9) for ^A a IJ and 
substitute the solution +4 A a IJ = + T a IJ back into (7.8), the additional equation of motion is 
automatically satisfied as a consequence of the Bianchi identity R[ a bc\d = 0. Hence, there are 
no "spurious" equations of motion. Moreover, since the self-dual and 3+1 Palatini actions 
differ by a term that is not a total divergence, the canonically conjugate variables for the two 
theories will disagree. As we shall see in the following section, it is this difference that will 
allow us to construct a Hamiltonian formulation of 3+1 gravity with a connection 1-form as 
the basic configuration variable. 

Finally, we conclude this subsection by showing the relationship between the self-dual 
and standard Einstein-Hilbert actions. To do this, note that since the equation of motion 
(7.9) for ' iA A a j J could be solved uniquely for ^A^ 3 in terms of the remaining basic variables 
4 e a , we can pull-back the self-dual action Ssd^^A) to the solution space +4 A a / J = + r a / J 
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and obtain a new action iS 5D ( 4 e). Doing this, we find 

_ 1 f ~abcd c A J i e J d i^L 

o JM 

where we expanded + R c d KL and used the Bianchi identity R[ a bc]d = to get the last line of 
(7.17). Thus, S_sd(^) is just 1/2 times the standard Einstein-Hilbert action SW( 4 e) viewed 
as a functional of a complex co-tetrad (See equation (6.5).) In fact, S_sd(S) = ^pO 5 )) 
where S_p{*e) is the pull-back of the 3+1 Palatini action. It was precisely to obtain this last 
equality that we defined the self-dual action (7.1) with an overall factor of 1/4 rather than 
1/8. 

7.2 Legendre transform 

To put the self-dual theory for complex 3+1 gravity in Hamiltonian form, we will basically 
proceed as we did in Section 6 for the 3+1 Palatini theory. However, since the spacetime 
metric g ab := ^ ^(fju is now complex, we can only assume that M is topologically Sxi? for 
some submanifold E and assume that there exists a real function t whose t = const surfaces 
foliate M. (We cannot assume that E is spacelike, since the signature of a complex metric 
is not defined.) We can still introduce a real flow vector field t a (satisfying t a (dt) a = 1) and 
a unit covariant normal n a to the t = const surfaces satisfying n a n a — — 1. (We are free to 
choose —1 for the normalization of n a since n a is allowed to be complex.) n a := g ab rib is the 
vector field associated to n a , and is related to t a by a complex lapse iV and complex shift 
iV a via t a = Nn a + iV a , with N a n a = 0. Finally, the induced metric q a b on E is given by 
q a b = g a b + n a n b . 

Following the same steps that we used in the previous section for the 3+1 Palatini theory, 
we find that (modulo a surface integral) the Lagrangian L SD of the self-dual theory is given 

b ^ , 

Lsd = / -NTr( + E a + E b + F ab ) + N a Tr( + E b + F ab ) 

Je (7.18) 
+ { + E a u )Cz + A a IJ + {^D a + E a jj)(^A ■ t) IJ . 

Here E a u denotes the self-dual part of E a jj := \tijKLX\ ahc ^ • Note that the transi- 
tion from the 3+1 Palatini Lagrangian to the self-dual Lagrangian can be made by simply 
replacing all of the real fields by their self-dual parts. The configuration variables of the 
theory are (^A -t) IJ , N, N a , + A a IJ , and + E a u . 

Now recall that for the real 3+1 Palatini theory, the configuration variable E a jj was not 
free to take on arbitrary values. From its definition in terms of the co-tetrad 4 e^, we saw that 

afe := e IJKL E a u E b KL = and Tr(E a E b ) > 0. (7.19) 
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The second condition followed from the fact that Tr(E a E b ) = 2q ab (= 2qq ab ), where q ab was 
the inverse of the induced positive-definite metric q a b on E. Taking the primary constraint 
(7.19) together with the 3+1 Palatini Lagrangian as the starting point for the Legendre 
transform, we found that the standard Dirac constraint analysis gave rise to additional 
constraints — one of which was 2nd class with respect to <fr ab = 0. By solving this 2nd class 
pair, the remaining (1st class) constraints became non-polynomial and we were forced back 
to the usual geometrodynamical description of real 3+1 gravity.^ 

Similarly, we must check to see if there are any primary constraints on the configuration 
variables of the self-dual theory. It turns out that although <p ab := e IJKL E a uE b kl — 
still follows from the definition of E a jj in terms of the complex co-tetrad it does not 
imply a constraint on the self-dual field + E a u. Equation (7.19) may be viewed, instead, 
as a constraint on the anti self-dual field ~E a jj. (Recall that for complex E a jj, ~E a jj is 
not necessarily the complex conjugate of + E a u). Thus, + E a u is free to take on arbitrary 
values, and the Legendre transform for the complex self-dual theory is actually fairly simple. 
By following the standard Dirac constraint analysis, we find that + E a ij is the momentum 
canonically conjugate to + A a IJ , while (^A ■ t) IJ , N, and iV a play the role of Lagrange 
multipliers. The complex phase space ( c Tsd, C ^sd) is coordinatized by the pair of complex 
fields ( + A a IJ , + E a u) and has the natural complex symplectic structure^ 

c n SD = I & + E a LJ A dt+4 a /J . (7.20) 

The Hamiltonian is given by 



H SD (+A, + E) = jf NTr( + E a + E b + F ab ) - N a Tr{+E b +F ab ) 

-{+D a + E a IJ ){+ A A-t) IJ . 



(7.21) 



As we shall see in the next subsection, this is just a sum of 1st class constraint functions 
associated with 

Tr( + E° + E b + F ab ) « 0, (7.22) 



It is fairly easy to see that all of the above statements — except for the non-holonomic constraint which 
would now say that Ti(E a E b ) be non-degenerate — apply to the complex 3+1 Palatini theory as well. <fi ab = 
is a primary constraint on the complex configuration variable E a u, and it must be included when performing 
the Legendre transform. The standard Dirac constraint analysis leads to a pair of 2nd class constraints which, 
when solved, gives back the usual geometrodynamical description of complex 3+1 gravity. 

25 Note that in terms of the Poisson bracket { , } defined by c Qsd, we have { + A a IJ (x), + E b KL(y)} = 
\5 b a 5(x,y){5^5 J I ] — %ckl MN $m$n)- The "extra" term on the right hand side is needed to make the Poisson 
bracket self-dual in the I J and KL pairs of indices. 
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Tr(+E b + F ab ) « 0, and (7.23) 
+D a + E a u « 0. (7.24) 

Note that all the constraints (and hence the evolution equations) are polynomial in the 
canonically conjugate pair ( + A a IJ , + E a u). This is a simplification that we found in the 
2+1 Palatini theory, but lost in the 3+1 Palatini theory when we solved the 2nd class 
constraints. In fact, since the constraint equations never involve the inverse of + E a u, the 
above Hamiltonian formulation is well-defined even if + E a u is non-invertible. Thus, we have 
a slight extension of complex general relativity. The self-dual theory makes sense even when 
the induced metric q ab = Tr( + E a + E b ) becomes degenerate. 

In order to make contact with the notation used in the literature (see, e.g., @), let us use 
the fact that the covariant normal n a to £ defines a unit internal vector rij via nj :— n a 
One can then show that 

+ b K u:=-~e K I j + iqf I n J] (7.25) 

is an isomorphism from the self-dual sub-Lie algebra of the complexified Lie algebra of 
5*0(3, 1) to the complexified tangent space of S. (Here e jkl '■= n^uxi, Qj '■= 5f + nin K , 
and rijn 1 — — 1.) It satisfies 

\ + b K u) ■= \{ + b K u - l -e u MN + b K MN ) = i + b K u, (7.26) 
[V, + b J } MN = e IJ K + b K MN , and (7.27) 
Tr(+6 J -V) := - + b\ IN + b JMN = -q IJ . (7.28) 

The inverse of + b K jj will be denoted by + bx IJ , and is obtained by simply raising and lowering 
the indices of + b K u with the internal metric rjjj. Since nx + b K u = 0, we will use a 3- 
dimensional abstract internal index i and write + b l u and + 6/ 3 in what follows. From property 
(7.27), it follows that + b l jj can actually be thought of as an isomorphism from the self-dual 
sub-Lie algebra of the complexified Lie algebra of SO(3, 1) to the complexified Lie algebra 
of SO{3). 

Given this isomorphism, we can now define a CCso (3)- valued connection 1-form A l a and 
a OC^^-valued vector density Ef via 

+ A a IJ =: A\ \ IJ and + E a u =: -iE? + b l u . (7.29) 

A straightforward calculation then shows thatQ 

\ 

26 To obtain equation (7.30), I assume that the fiducial derivative operator d a has been extended to act on 
C£so(3)- m dices in such a way that d a + bi IJ = 0. 



IJ = (2d [a Ai + e l jk A{A k b ) V J =: F\ h \ IJ and (7.30) 
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Tr( + E a + E b ) = E a i E u = <f 6 '. (7.31) 

Thus, Ef is a complex (densitized) triad and F % ah is the Lie algebra-valued curvature tensor 
of the generalized derivative operator T> a defined by T> a v % := d a v l + t l jkA J a v k '. 
In terms of A\ and Ef, the complex symplectic structure c flsD becomes 

C Q SD = -i [ dEf A cL4* , (7.32) 

so — i-E 1 ? is the momentum canonically conjugate to A\. The Hamiltonian (7.21) can be 
written as 

H SD (A, E) = jf ^ Ne ijk EfE b F abk - iiV°^i^ + i(P«£?)( 4 .A • t)\ (7.33) 

while the constraint equations (7.22)-(7.24) can be written as 

e^ k EfE b F abk » 0, (7.34) 
^i&«0, and (7.35) 
P a ^ a » 0. (7.36) 

We will take the constraint equations in this form when we analyze the Poisson bracket 
algebra of the corresponding constraint functions in the following section. 

So far, all of the discussion in this section has dealt with complex 3+1 gravity. In order 
to recover the real theory, we must now impose reality conditions on the complex phase space 
variables (A\, Ef) to select a real section of ( c T S d, c ^sd)- To do this, recall that in terms of 
the standard geometrodynamical variables (q a b,P ab ), one recovers real general relativity from 
the complex theory by requiring that q ab and p ab both be real. Since equation (7.31) tells us 
that E?E bl = q ab (= qq ab ), the condition that q ab be real can be conveniently expressed in 
terms of Ef as 

EfE u be real. (7.37) 

Since we will want to ensure that this reality condition be preserved under the dynamical 
evolution generated by the Hamiltonian, we must also demand that 

(EfE bi y be real. (7.38) 

Since in a 4-dimensional solution of the field equations p ab is effectively the time derivative 
of q ab , requirement (7.38) is equivalent to the condition that p ab be real. In addition, since 
the Hamiltonian of the theory is just a sum of the constraints (7.34)-(7.36) (all of which are 
polynomial in the canonically conjugate variables), the reality conditions (7.37) and (7.38) 
are also polynomial in (A l a ,Ef). 
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Finally, to conclude this subsection, I should point out that the self-dual action (7.1) 
viewed as a functional of a self-dual connection 1-form +i A aI J and a real co-tetrad 4 e^ does 
not yield the new variables for real 3+1 gravity when one performs a 3+1 decomposition. 
The definition of the configuration variable + E a u in terms of a real co-tetrad 4 e^ gives rise 
to a primary constraint. Although the non-holonomic constraint can be expressed in terms 
of + E a LJ as 

Tr(+E a + E b ) > 0, (7.39) 

the holonomic constraint <p ab = cannot be expressed solely in terms of + E a u. For real E a u 
we have that ~E a u equals the complex conjugate of + E a u, so 



<T = e IJKL ( + E a jj + +E» IJ ){ + E b KL + + E b KL ) = 0. (7.40) 

But by writing 4> ab = in this way, we have destroyed the possibility of completing the 
standard Dirac constraint analysis. For nowhere in the analysis have we been told how to 
take Poisson brackets of the complex conjugate fields. The Legendre transform of the self- 
dual Lagrangian for real 3+1 gravity breaks down when we try to incorporate the primary 
constraints into the analysis. 

7.3 Constraint algebra 

Given constraint equations (7.34)-(7.36) for the complex self-dual theory, we would now 
like to verify the claim that their associated constraint functions form a 1st class set. To 
do this, let v l (which takes values in CCso(3)), N, and N a be arbitrary complex- valued test 
fields on E and define 

C(N):=y^ k E?E b E abk , (7.41) 
C\N) := —iJ^N a E b F i ah , and (7.42) 
G(v) := -i [ v\V a Et). (7.43) 



E 

These will be called the scalar, vector, and Gauss constraint functions. As the names and 
notation suggest, these constraint functions will play a similar role to the constraint functions 
defined in subsection 5.3. Many of the calculations and results found there will apply here 
as well. 

As usual, it is fairly easy to show that the Gauss constraint functions generate the 
standard gauge transformations of the connection 1-form and rotation of internal indices. 
Since 

6G{V) ^V a v* and ^ = E% (:= -if^E*) , (7.44) 



6E? SA 



a 
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it follows that A\ i-> A\ - eV a v l + 0(e 2 ) and Ef i-> £7? - e{v, E a }i + 0(e 2 ). It also follows 
that 

{G(v),G(w)} = G([v,w}), (7.45) 

where := e % jkV^w k is the Lie bracket of v % and w l . Thus, the mapping v i— > is a 

representation of the Lie algebra CCso(3)- Furthermore, given its geometrical interpretation 
as the generator of internal rotations, we have 

{G(v),C(N)} = and (7.46) 
{G(v),C'(N)} = 0, (7.47) 

as well. 

Since it is possible to show that the vector constraint function does not by itself have any 
direct geometrical interpretation (see, e.g., |34]| ) we will define a new constraint function, 
C(N), by taking a linear combination of the vector and Gauss constraints. We define 

C(N) := C'{N) - G(N), (7.48) 

where N l := N a A l a . We will call C(N) the diffeomorphism constraint function since the 
motion it generates on phase space corresponds to the 1-parameter family of diffeomorphisms 
on S associated with the vector field N a . To see this, we can write 

C{N) : = C\N) - G(N) 

= -i( N a E^F^ b + i [ N\V a E«) 

= -i / N^E F^-AiV^) 

= -i jf EtC^Ai 

where the Lie derivative with respect to iV a treats fields having only internal indices as 
scalars. To obtain the last line of (7.49), we ignored a surface integral (which would vanish 
anyways for iV a satisfying the appropriate boundary conditions). By inspection, it follows 
that A % a I— >■ A\ + eC^A l a + 0(e 2 ), etc. Using this geometric interpretation of C(N), it follows 
that 

{C(N),G(v)} = G(C$v), (7.50) 
{C(N),C(M)} = C(C^M), and (7.51) 
{C{N),C{M)} = C{[N,M\). (7.52) 
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We are left to evaluate the Poisson bracket {C(N),C(M)} of two scalar constraints. 
Using 

6C(N) 6C(N) .. _ ~. 

~ Ne^ k E]F ahh and — ^ = e^ k V b (NE^E b ), (7.53) 



(1 



it follows that 



r SC(N) 6C(M) 

{C(N) J C(M)}= / -¥r- — (N<-+ M) 

w V ~ ;J 7e 54^ 5(-iE a ) ~ 



/ ier^NE^Me^E^ - (N <- M) (7-54) 
ie^er n KE c n E b JMd c N - Nd c M)F abk . 



If we now use the fact that 

(which is a property of the structure constants of SO (3)), we get 

{C(N),C(M)} = C\K) (= C(K) + G{K)) , (7.56) 

where K a := q ab (Nd b M — Md b N) and g afc = E^E b \ Thus, the constraint functions are 
closed under Poisson bracket — i.e., they form a 1st class set. Note, however, that since the 
vector field K a depends on the phase space variable Ef, the Poisson bracket (7.56) involves 
structure functions. The constraint functions do not form a Lie algebra. 



8. 3+1 matter couplings 

In this section, we will couple various matter fields to 3+1 gravity. We will repeat much of 
what we did in Section 5, but this time in the context of the 3+1 theory, and for a Yang-Mills 
field instead of a massless scalar field. In subsections 8.1 and 8.2, we couple a cosmological 
constant A and a Yang-Mills field to complex 3+1 gravity using an action principle and the 
self-dual action as our starting point. We shall show that the inclusion of these matter fields 
does not destroy the polynomial nature of the constraint equations. This is the main result. 
(As usual, reality conditions should be included to recover the real theory.) As I mentioned 
for the 2+1 theory, it is possible to couple other fundamental matter fields (e.g., scalar and 
Dirac fields) to 3+1 gravity in a similar fashion and obtain the same basic results. For a 



more detailed discussion of this and related issues, interested readers should see, e.g., 
8.1 Self-dual theory coupled to a cosmological constant 
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To couple a cosmological constant A to complex 3+1 gravity via the self-dual action, we 
will start with the action 

S^A) := \ f M V abcd e IJK L V a V b (V L - A %^ e j). ( 8 .1) 

Here ^F abI J = 2d[ a ~ +4 A b ] I J + \^A a , ^4 fe ]/ J is the internal curvature tensor of the self-dual 
generalized derivative operator JrA D a defined by the self-dual connection 1-form ^Aaf 1 , and 
^ is a complex co-tetrad which defines a spacetime metric g ab via g ab :— ^e a ^{vu- Note 
that Sa^, is just a sum of the self-dual action 

S SD ( 4 e, «A) :=-( V abcd euKL H V b ^F cd KL (8.2) 

4 JM 

and a term proportional to the volume of the spacetime. In fact, 

77 / ?7 e/XKL e a e fe e c e d — A V~9, (s.dj 

4! Jm jM 

where g is the determinant of the covariant metric g ab . 

To show that (8.1) reproduces the standard equation of motion, 

G ab + Ag ab = 0, (8.4) 

for gravity coupled to the cosmological constant A, we will first vary (8.1) with respect to 
the self-dual connection 1-form +4 A</ J '. Since the second term (8.3) is independent of ^Aj 3 , 
we get 

(( 4 e) + ( 4 4° 4 e 6 J)) = 0, (8.5) 

which is exactly the equation of motion we obtained in Section 7 for the vacuum case. Thus, 
just as we saw in subsection 7.1, ^A^ 3 = + T a j J where + T a i J is the self-dual part of the 
internal Christoffel symbol of V a (the unique, torsion-free generalized derivative operator 
compatible with the co-tetrad.) Since (8.5) can be solved uniquely for ^A^ 3 in terms of 
the remaining basic variables *e a , we can pull-back S^{*e, ^A) to the solution space +4 A a j J = 
+ T aI J . We obtain a new action 

£A( 4 e) = \ J m V abcd e IJKL V a V b (+R cd KL - A %f 4 e^), (8.6) 

where + R cd KL is the self-dual part of the internal curvature tensor defined by V a . Then by 
using the Bianchi identity R[ a bc]d = for the first term and (8.3) for the second, we get 

S^e) = \ I J=j(R-2A). (8-7) 
2 jm 
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As mentioned in subsection 5.1, this is (up to an overall factor of 1/2) the action that one 
uses to obtain (8.4) starting from an action principle. This is the desired result. 

To put this theory in Hamiltonian form, we proceed as in subsection 7.2. Recall that 
(modulo a surface integral) the Lagrangian L SD of the self-dual theory is given by 

Lsd = I -NTr( + E a + E b + F ab ) + N a Tr( + E b + F ah ) 

Je (8.8) 
+ ( + E a u )£ t +A a IJ + {+D a + E a u )^A ■ t) IJ , 

where + E a u denotes the self-dual part of E a u := \eijKLV abc 4e b ■ By using the isomor- 
phism between the self-dual sub-Lie algebra of the complexified Lie algebra of 5*0(3, 1) and 
the complexified Lie algebra of S0(3), we can rewrite (8.8) as 

Lsd = jf — ^ Ne ijk E^EjF abk + iN a E b F^ - iE?£ t A\ - i(V a E%)?A ■ t)\ (8.9) 

where Ef is a complex (densitized) triad (i.e., E£E bl = q ab (= qq ab )) and A l a is a connection 
1-form on X that takes values in the complexified Lie algebra of S0(3). 
By using the decomposition y/^g = N^/q dt together with the fact that 

^ bc e^ k EtE b E c k = q, (8.10) 



one can similarly show that 

A I" ^,„r,„ r l ATA1AISAT A 

4\Jm' 

Thus, the Lagrangian L A for 3+1 gravity coupled to the cosmological constant A via the 
self-dual action is given by 



J V abcd euKL V a V b 4 ef % L d = A J dt j N^ k E?E b E c k . (8.11) 



= I -N{\t ijk E?E]F abk + ^ bc e« k E«E]Ei) 

+ iN a E b F l ab - iE?£f A\ - t(V a E?)( A A ■ t)\ 

The configuration variables of the theory are ( A A ■ t) 1 , N, N a , A l a , and Ef. 

By following the standard Dirac constraint analysis, we find (as in the vacuum case) that 
—iE£ is the momentum canonically conjugate to A l a while ( 4 A -t) 1 , N, and N a play the role 
of Lagrange multipliers. The complex phase space and complex symplectic structure are the 
same as those found for the self-dual theory with A = 0, while the Hamiltonian is given by 

Ha(A E) = jf N^ k EtE b F abk + ^e^E^Et) 

-zN a E b F* b + z(V a E?)CA-ty. 
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We shall see that this is just a sum of 1st class constraint functions associated with 

^ k E?E b F abk + ^ hc ^ k E<tE)El « 0, (8.14) 

W«0. and (8.15) 
£> a £f«0. (8.16) 

These are the constraint equations associated with the Lagrange multipliers N, N a , and 
( 4 A ■ t) 1 , respectively. Note that they are polynomial in the canonically conjugate variables 
(A l a ,E£) even when A ^ 0. In fact, only constraint equation (8.14) differs from its A = 
counterpart. 

To conclude this subsection, we will verify the claim that the constraint functions associ- 
ated with (8.14)-(8.16) form a 1st class set. Since the Gauss and diffeomorphism constraint 
functions associated with (8.16) and (8.15) will be the same as in subsection 7.3, we need 
only concentrate on the scalar constraint function 

C(N) := j^N^ k E?E b F abk + ^e^E^E*). (8.17) 



Since 



SC(N) A 

-j^ = N^ k E b F abk + -^ k E)E c k ) and (8.18a) 
SC(N) 



SAi 



= ei k V b {NE«El), (8.186) 



it follows that 



r 5C(N) 5C(M) 
{C(N), C(M)} = / — ^ ^— - (N <-> M) 

A 



ter n V c {NE a m El)M^ k E b 3 F abk + -^ k E b E d k ) - (N <- M) 
= J^ k er n {Md c N - Nd c M)(E a m E c n E b F abk + ±qe mjk E c n ). 

(8.19) 

If we again use the fact that the structure constants of 5*0(3) satisfy 

e ijk e mn = ( 5 jm 5 kn _ gjngkm^ (g^Q) 

we get 

{C(N),C(M)} = C'(K) (=C(K) + G(K)), (8.21) 

where K a := q ab (Nd b M — Md b N) and q ab = E£E bl as before. Thus, the constraint functions 
are closed under Poisson bracket — i.e., they form a 1st class set. The Poisson bracket algebra 
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of the constraint functions is exactly the same as it was for the A = case. In particular, 
since the vector field K a depends on the phase space variable Ef, the constraint functions 
again do not form a Lie algebra. 

8.2 Self-dual theory coupled to a Yang-Mills field 

To couple a Yang-Mills field (with gauge group G) to complex 3+1 gravity via the self- 
dual action, we will start with the total action 

S T {% % 4 A) := S SD ( 4 e, «A) + \s YM {% 4 A), (8.22) 

where SsdO 2 , ~^A) is the self-dual action (8.2) and SymC^, 4 A) is the usual Yang-Mills action 

SymCc 4 A) := - / Tr(v^ g ac g bd *F ab *F cd ). (8.23) 

JM 

Here Sym{^i 4 A) is to be viewed as a functional of a co-tetrad ^ and a connection 1-form 4 A a 
which takes values in the Lie algebra of the gauge group G.Q Tr denotes the trace operation 
in some representation of the Yang-Mills Lie algebra, and ^ ab = 2<9[ a 4 A fe ] + [ 4 A a , 4 Ab] is the 
(internal) curvature tensor of the generalized derivative operator ^Da defined by 4 A a . The 
additional factor of 1/2 is needed in front of S , yjv/( 4 e, 4 A) so that the above definition of 
the total action will be consistent with the definition of Ssd{ 4 £, a A). The Yang-Mills action 
depends on the co-tetrad ^ through its dependence on y/^g and g ab , but is independent of 
the self-dual connection 1-form 44 A a / J . As mentioned in Section 5, out of all the fundamental 
matter couplings, only the action for the Dirac field would depend on 44 A a / J . 

To show that (8.22) reproduces the standard Yang-Mills coupled to gravity equations of 
motion 

*D b (V=g 4 F a6 ) = and G ad = 87cT ad {YM), (8.24) 

where 

T ab (YM) := -j-Trfa* *F bc - ^ ab 4 F C(i 4 F cd ) (8.25) 

is the stress-energy tensor of the Yang-Mills field, we proceed as we did in the previous 
subsection. Since S'yjv ! f( 4 e, 4 A) is independent of ~^A a IJ , the variation of (8.22) with respect 
to A a implies 

*V h (( 4 e) + ( 4 ei a VJ)) = 0. (8.26) 

27 Yang-Mills fields will be denoted by bold face stem letters and their (internal) Lie algebra indices will 
be suppressed. Throughout, we will assume that we have a representation of the Yang-Mills Lie algebra Cg 
by linear operators (on some vector space V) with the trace operation Tr playing the role of an invariant, 
non-degenerate bilinear form k. 
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As before, this tells us that ^Aai 3 = + T aI J . Recalling that the Bianchi identity R[ a bc]d = 
implies that the pull-back of SsdC^, ~^A) to the solution space ^A^ 3 = + T aI J is just 1/2 
times the standard Einstein-Hilbert action Seh^) f° r complex 3+1 gravity, we obtain 

S T ( 4 e, 4 A) = ^(S EH ( 4 e) + S yM ( 4 e, 4 A)). (8.27) 

This is (up to an overall factor of 1/2) the usual total action that one uses to couple a Yang- 
Mills field to gravity. If we now vary 5t(%, 4 A) with respect to 4 A a and 4 e a , and contract the 
second equation with "fef, we recover (8.24). Note that to write the first equation in (8.24), 
we had to consider a torsion-free extension of ^Da to spacetime tensor fields. But since the 
left hand side is the divergence of a skew spacetime tensor density of weight +1 on M, it is 
independent of this choice. 

To put this theory in Hamiltonian form, we need only decompose the Yang-Mills action 
Sym^, 4 A) since the self-dual Lagrangian L SD is given by (8.9). Using g ab = q ab — n a n b and 
•>/— p = N^/q dt it follows that 

S Y m(WA) = I dt I Tr{ - Nq- 1 r c f d F ab F cd + 2r b N' 1 q- 1 x 

J Js L (8.28) 
x (£ r A a - D a ( 4 A • t) + iV c F ac )(£ f A 6 - D 6 ( 4 A • t) + N d F bd )}, 

where <f fe := qq ab (= EfE"), ( 4 A • t) := t a 4 A a , and A a := q b 4 A b . Here F ab := q c a q d4 F cd is 
the curvature tensor of the generalized derivative operator D a (:= q b ^D&) on £ associated 
with A a . If we now define the "magnetic field" of A a to be B ab := 2F ab (= 2q c a q d ^cd), we 
see that the Yang-Mills Lagrangian L Y m is given by 

L YM = [ Tr{ - l -Nq- l r c f d B ab B cd + 2? b N- 1 q- 1 x 

, , (8-29) 

x (£ t - A a - D a ( 4 A • t) + -iV c B ac )(£ r A 6 - D 6 ( 4 A • t) + -N d B bd )}. 

The total Lagrangian Lt is the sum = L$d + \Lym and is to be viewed as a functional of 
the configuration variables ( 4 A-t), ( 4 A-ty, N, N a , A a , Ef, A a and their first time derivatives. 
Following the standard Dirac constraint analysis, we find that 

E> := ^ = 2? b N- 1 q-\jC r A b - D b ( 4 A • t) + l -N d B bd ) (8.30) 

is the momentum (or "electric field") canonically conjugate to A . Since this equation can 
be inverted to give 

C r A a = ^q ab NE b + D a ( 4 A • t) - ^N c B ac , (8.31) 
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it does not define a constraint. On the other hand, —iEf is constrained to be the momentum 
canonically conjugate to A\, while ( 4 A • t), ( 4 A ■ t)\ N, and N a play the role of Lagrange 
multipliers. The resulting complex total phase space ( C T T , c f2 T ) is coordinatized by the 
pairs of fields (A a) E£) and (A a , E a ) with symplectic structure 

c Vt T = -i&E? A &AI + Tr(dtE a A dA ). (8.32) 

The Hamiltonian is given by 

H T (A, E, A, E) = N^ k E?E]F abk + ^ W d Tr(B ab B cd + E ab E cd )) 
+ N a [ - iE\F l ah + Tr(E 6 F a6 )) + i(V a E°)( AA ■ t) 1 - Tr(( 4 A • t)D a E a ), 



(8.33) 



where E a b := 7^fe c E c is the dual to the Yang-Mills "electric field" E a . We shall see that this 
is just a sum of 1st class constraint functions associated with 

1 e ijkE-E b jFahk + ± q -^ c f d Tr(B ab B cd + E ab E cd ) » 0, (8.34) 
- iE\Fl h + Tr(E b F a6 ) » 0, (8.35) 
V a E1 sa 0, and D a E a w 0. (8.36) 

These are the constraint equations associated with the Lagrange multipliers N, N a , ^A-t) 1 , 
and ( 4 A • t) , respectively. 

Note that by inspection (8.35) and (8.36) are polynomial in the canonically conjugate 
variables. However, constraint equation (8.34) fails to be polynomial due to the presence of 
the non-polynomial multiplicative factor q~ x . But since q = ^'fl^e t i k E£EjE%. is polynomial 
in Ef, we can multiply (8.34) by q and restore the polynomial nature of all the constraints. 
Thus, to couple a Yang-Mills field to 3+1 gravity via the self-dual action, we are led to a 
scalar constraint with density weight +4. This implies that the associated constraint function 
will be labeled by a test field (i.e., lapse function) having density weight —3. 

To verify the claim that the constraint functions associated with (8.34)-(8.36) form a 
1st class set, let v l and v (which take values in complexified Lie algebra of SO (3) and the 
representation of the Lie algebra of the Yang-Mills gauge group G), N, and N a be arbitrary 
complex- valued test fields on E. Then define 

C(N) := jT N{^ k E1E]F abk + ^~ W d Tr(B ab B cd + E ab E cd )) , (8.37) 

C\N) := N a (-iE^L + Tr(E 6 F afc )), and (8.38) 

G(v, v) := j s Tr(vD a E a ) - iv^Ef), (8.39) 
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to be the scalar, vector, and Gauss constraint functions. 

As usual, it is fairly easy to show that the Gauss constraint functions generate the 
standard gauge transformations of the connection 1-forms and rotations of internal indices. 
Using this information, we find that 

{G(v, v), G(w, w)} = G([v, w], [v, w}), (8.40) 
{G(v,v),C(N)} = 0, and (8.41) 
{G(v,v),C'(N)} = 0, (8.42) 

where [v,w] and w] 1 are the Lie brackets in Cg and CCso{z)- Thus, the subset of Gauss 
constraint functions form a Lie algebra with respect to Poisson bracket. In fact, the mapping 
(v, v) I— > G(v,v) is a representation of the direct sum Lie algebra Cg © CC S o(3)- 

Again, the the vector constraint function will not have any direct geometrical inter- 
pretation, so we define the diffeomorphism constraint function C(N) by taking a linear 
combination of the vector and Gauss law constraints. Setting 

C(N) := C\N) - G(N, N), (8.43) 

where N := N a A a and N l := N a A l a , we can show that 

C(N) = [ -iE^nK + TME^AJ, (8.44) 

where the Lie derivative with respect to A^ a treats fields having only internal indices as 
scalars. By inspection, A l a i— > A l a + eC^A l a + 0(e 2 ), etc., so the motion on phase space 
generated by C (N) corresponds to the 1-parameter family of diffeomorphisms on E associated 
with A^ a . From this geometric interpretation of C(N), it follows that 

{C(N), G(v, „)} = G(C M v, Cjjv), (8.45) 
{C(N),C(M)} = C(C^M), and (8.46) 
{C(N), C(M)} = C([N, M\). (8.47) 

Finally, we are left to evaluate the Poisson bracket {C(N),C(M)} of two scalar constraint 
functions. After a fairly long but straightforward calculation that uses the fact that the 
structure constants of 5*0(3) satisfy 

e ijk e mn = ( S jm 5 kn _ gjngkm^ ( 8>48 ) 

one can show that 

{C(N),C(M)} = C'(K) (=C(K) + G(K,K)), (8.49) 
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where K a := q ab (NdbM — MdbN) and q ab = EfE bl . Thus, the constraint functions are again 
closed under Poisson bracket — i.e., they form a 1st class set. Just as we saw in subsection 
8.1 for the cosmological constant A, the Poisson bracket algebra of the constraint functions 
for complex 3+1 gravity coupled to a Yang-Mills field via the self-dual action is exactly the 
same as it was for the vacuum case. 

9. General relativity without-the-metric 

To conclude this review, we will describe a theory of 3+1 gravity without a metric. 
This will complete the transition from geometrodynamics to connection dynamics in 3+1 
dimensions. Although we saw in Section 7 that the Hamiltonian formulation of the self-dual 
theory for complex 3+1 gravity could be described in terms of a connection 1-form A\ and its 
canonically conjugate momentum (or "electric field") Ef, the action for the self-dual theory 
depended on both a self-dual connection 1-form +4 A ( / J and a complex co-tetrad ^e l a . Since 
the co-tetrad defines a spacetime metric g a b via g a b '■= 4 e« ^(rju, the self-dual action had an 
implicit dependence on g a b. The purpose of this section is to show that (modulo an important 
degeneracy) complex 3+1 gravity can be described by an action which does not depend on 
a spacetime metric in any way whatsoever. We shall see in subsection 9.1 that this action 
depends only on a connection 1-form A A l a (which takes values in the complexified Lie algebra 
of 5*0(3)) and a scalar density $ of weight -1 on M. Hence we obtain a pure spin- connection 
formulation of gravity. We shall also see how this pure spin-connection action is related to 
the self-dual action in the non-degenerate case. 

In subsection 9.2, we will analyze the constraint equations for this theory. Since we will 
have shown in subsection 9.1 that the self-dual action and the pure spin-connection action 
are equivalent when the self-dual part of the Weyl tensor is non-degenerate, the constraint 
equations of this theory are the same as the the constraint equations for the self-dual theory 
found in subsection 7.2. However, we will now be able to write down the most general 
solution to the four diffeomorphism constraint equations (the scalar and vector constraints 
of the self-dual theory) when the "magnetic field" Bf associated with the connection 1-form 
A l a is non-degenerate. This is a new result for the Hamiltonian formulation of the 3+1 theory 
that was made manifest by working in the pure spin-connection formalism. 

I should emphasize here that all of the results in this section are taken from previous 
work of Capovilla, Dell, Jacobson, Mason, and Plebanski. I am not adding anything new in 
this section; rather, I am reporting their results to bring the discussion of geometrodynamics 
versus connection dynamics for 3+1 gravity to it logical conclusion. Readers interested in a 
more detailed discussion of the general relativity without-the-metric theory (including matter 
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couplings and an extension of this theory to a class of generally covariant gauge theories) 



should see |T0| , |1T| , |12| , |13| , |14l and references mentioned therein. In addition, Peldan has 
recently provided a similar pure spin-connection formulation of 2+1 gravity. Interested 
readers should see [|1^ . 

9.1 A pure spin-connection formulation of 3+1 gravity 

The pure spin- connection action for complex 3+1 gravity is defined to be 

S(& A A) :=\l %(rj. 4 F* A 4 F% • 4 F k A 4 F<>^, (9.1) 

o JM 

where is a connection 1-form which takes values in the complexified Lie algebra of 5*0(3), 
$ is a scalar density of weight -1 on M, and hijki and (fj- A F % A 4 F J ) are shorthand notations 
for 

%fc; := (^ifc^ji + 8u$jk - and (9.2) 

(fj- <PA 4 ^):=^ 4 F^. (9.3) 

As usual, 4 F ( J 6 = 29[ a + [ A OI is the Lie algebra-valued curvature tensor of the gener- 
alized derivative operator *D a defined by 

*V a rf:=d a rf + \ A A ai v\\ (9.4) 

where [ 4 v4. a , u]' := e l j k 4 A 3 a v k denotes the Lie bracket of 4 A l a and v l in CCso(3)- Although *D a 
defined by (9.4) knows how to act only on internal indices, we will often find it convenient 
to consider a torsion-free extension of ^Va to spacetime tensor fields. All results and all 
calculations will be independent of this choice. 

To show that the pure spin-connection action reproduces the standard results of complex 
3+1 gravity, one could vary (9.1) with respect to $ and A A l a and analyze the resulting Euler- 
Lagrange equations of motion. Instead, we will start with the self-dual action 

S SD ( 4 c, «A) := \ I V abcd e IJK L A e{ *e J b V L (9.5) 

for complex 3+1 gravity and show that (modulo an important degeneracy) the self-dual 
action (9.5) is actually equivalent to (9.1). Basically, we will eliminate from (9.5) the field 
variables which pertain to the spacetime metric by solving their associated Euler-Lagrange 
equations of motion. This will require that a certain symmetric trace-free tensor ipij be 
invertible as a 3 x 3 matrix. By substituting these solutions back into the original action 
(9.5), we will eventually obtain (9.1). We should point out that, in a solution, ipij corresponds 
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to the self-dual part of the Weyl tensor associated with the connection 1-form A A l a . Thus, the 
equivalence between the two actions breaks down whenever the self-dual part of the Weyl 
tensor is degenerate. Note also that the pure spin-connection action describes complex 3+1 
gravity. To recover the real theory, one would have to impose reality conditions similar to 
those used in Section 7 for the self-dual theory. For a detailed discussion of ipij and the 
reality conditions see, e.g., [O. 

Since the self-dual action (9.5) depends on both a self-dual connection 1-form ^Aj -1 and a 
complex co-tetrad it has an implicit dependence on the spacetime metric g ab := *e l a ^rjij. 
Thus, it should come as no surprise that the first step in obtaining a metric-independent 
action for 3+1 gravity involves the elimination of ^ from (9.5). To do this, let us define 

Sab/j := eijKL 4 ef A e b (9.6) 
and + S a b/j to be its self-dual part.0 Then we can write the self-dual action as 

5 SD ( 4 e, «A) = \ I v abcd ^a b u ^F cd IJ , (9.7) 
4 Jm 

where we have used the fact that T> abIJ ^F cd IJ = + T> a bu ^F crf /J . To simplify the notation 
somewhat, let us recall that the self-dual sub-Lie algebra of the complexified Lie algebra 
of 5*0(3, 1) is isomorphic to the complexified Lie algebra of 5*0(3). Using the isomorphism 
described in Section 7, we can define a O/^o^-valued connection 1-form 4 A a and a CC* so( ^ 3 y 
valued 2-form T, abi via 

V=: 4 4V J and ^ abIJ =: E abi (9.8) 

Then 

S SD { A eM) = - [ fj abcd E abi 4 F* d , (9.9) 

4 J M 

where A F % ab = 2d[ a t4i + e % ^ A A{ 4 A b is the Lie algebra-valued curvature tensor of the gener- 
alized derivative operator *D a defined by 4 A a . It is related to ^ A F ab IJ via ^F^ 1,1 = 4 F^ b + bi IJ . 

Although the right hand side of (9.9) involves just S afc j and 4 A l a , the action is still a 
functional of 4 A l a and since T* a bi depends on 4 e a through equation (9.6). To eliminate ^ 
from the action, we must use the result (see, e.g., |12| ) that (9.6) holds for some 4 e a if and 
only if the trace-free part of E* A S J equals zero — i.e., £ a &/j = €jjkl Ae b f° r some if 
and only if 

rj^&A ~ l^a^dk) = 0. (9.10) 

28 Recall that the self-dual part of XWj is defined by + S a fc/j := h(^abU ~ h e U KL ^abKL)- 
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Thus, the self-dual action can be viewed as a functional of instead of provided we 
include in the action a term which gives back (9.10) as one of its Euler-Lagrange equations 
of motion. More precisely, let us define 

S, 4 A) :=- ( r?**p M 4 F* cd - ^KA), (9-H) 

4 JM I 

where ip^ is a symmetric trace-free tensor which will play the role of a Lagrange multiplier 
of the theory. Then the variation of S(ip, £, 4 -A) with respect to ipij will yield (9.10). Solving 
this equation and pulling-back the action (9.11) to this solution space gives back (9.9). 

But instead of varying S(ip, S, 4 A) with respect to ipy, let us vary this action with respect 
to S a fej and solve the resulting Euler-Lagrange equation of motion for E a &j in terms of ipij 
and 4 A a . Varying (9.11) with respect to T, a bi, we find 



i 



^-#^ = 0, (9.12) 

where ip % i := 5 S^ip^. This equation can be solved for £ a fo in terms of the remaining field 
variables provided the inverse {ip~ l )ij of ipy exists. Assuming that it does, we get 



J abi 



(r% 'Kb- (9.i3) 



If we now pull-back (9.11) to the solution space defined by (9.13), the resulting action 
becomes 

S(^A) = \l V abc \r% A Fl b 4 Fl d . (9.14) 

o JM 

This is to be viewed as a functional of only the symmetric trace-free tensor ip^ and the 
connection 1-form 4 A a . 

We are almost finished. What remains to be shown is that ip^ can be eliminated from 
the action (9.14) in lieu of a scalar density $ of weight -1 on M. To do this, let us write the 
action in matrix notation and introduce another Lagrange multiplier fx to guarantee that ipij 
is trace-free.0 Then (9.14) can be written as 

Sifi^M) = \ I Tr(-0 _1 M) + A*13r0, (9.15) 

8 JM 

where ipij is now assumed to be only symmetric (and invertible) and M 4 - 7 is defined by 

M ij := fj abcd *Fl h *Fi d . (9.16) 
Varying S(J1, ip, A A) with respect to ipij, we find 

—ip^Mip' 1 +/1I — 0. (9.17) 

29 By introducing Jl, we can consider arbitrary symmetric variations of ipij rather than symmetric and 
trace-free variations. 
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By multiplying on the left and right by ip, we see that (9.17) is equivalent to 

M = Jii) 2 . (9.18) 

This equation can be solved for ipij in terms of and Jl provided the square-root of 
exists. Then 

i> = ji- l/2 M 1 / 2 , (9.19) 
so the action (9.15) pulled-back to this solution space equals 

S(ji, *A) = ^J Ji l ' 2 TiM 1 ' 2 . (9.20) 

The variation of S(J1, A A) with respect to Jl now implies that TrM 1 / 2 = 0. From (9.19) we 
see that this is nothing more than the tracelessness of ipij . 

In order to write the action in its final form (9.1), recall that the characteristic equation 
obeyed by any 3x3 matrix is 

B 3 - (TrB)B 2 + - ((TrB) 2 - TrB 2 ) B - (detB)J = 0. (9.21) 

Multiplying by B and setting B 2 = M (i.e., B = M 1/2 ), we get 

(detB)B = M 2 — -(TrM)M, (9.22) 

(Here we have used the fact that TrB (= TrM 1 / 2 ) = 0.) Using (detB) 2 = detM and 
assuming invertibility of B (so that detB ^ 0), we can write this last equation as 

B = (detM)~ 1/2 (M 2 - -(TrM)M). (9.23) 

By substituting this expression for B (= M 1//2 ) back into (9.20), we find 

S(fi,*A) = 7 / /I 1/2 (detM)- 1/2 Tr(M 2 - -(TrM)M). (9.24) 
4 Jm 2 

Finally, if we define 

$ = /I 1/2 (detM)- 1/2 (9.25) 

(which is a scalar density of weight -1 on M) and use the definition (9.16) of in terms 
of 4 F* 6 , we see that 

S{%M) = \l m- ^ A 4 F%- 4 F k A 'F l )h l3kl (9.26) 

o JM 

when viewed as a functional of $ and 4 A l a . Note that hijki and (rj ■ A F % A 4 F J ) are given as 
before by equations (9.2) and (9.3). This is the desired result. 
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9.2 Solution of the diffeomorphism constraints 

Given that the self-dual and pure spin-connection actions are equivalent when the self- 
dual part of the Weyl tensor is non-degenerate, it follows that the constraint equations of 
the theory can be written as 

e ijk^b Fgbk w 0> (a27) 

£?i&»0, and (9.28) 
V a E1 « 0. (9.29) 

These are just the constraint equations that we found in subsection 7.2 when we put the 
self-dual theory in Hamiltonian form. As before, the canonically conjugate variables consist 
of a pair of complex fields (A l a , Ef), where A l a is the pull-back of the connection 1-form A l a 
to the submanifold S and Ef is a complex (densitized) triad which may or may not define 
an invertible induced metric q ab := EfE bl . However, by working in the pure spin-connection 
formalism, we will obtain a new result. We will be able to write down the most general 
solution to the four diffeomorphism constraints (9.27)-(9.28) when the "magnetic field" Bf 
associated with A l a is non-degenerate. 

To see this, recall that in the self-dual theory 

+ E a LJ =: -iE? (9.30) 

where E a jj was the self-dual part of 

E a u := l -e U KLV abc 4 ef 4 e c L . (9.31) 

Note that in terms of T> abIJ defined by (9.6), we have E a u = ^fj abc T, bcIJ , so that 

-iEf = \v abc ^bd. (9.32) 
If we now use the result that an invertible symmetric trace-free tensor ipij implies 

Z abl = A Fl h) (9.33) 

it follows that 

Ef = i{r%B a \ (9.34) 

where B m := \r\ ahc 4 F bc (= |?7 abc F b * c ) is the "magnetic field" of A l a . We will now show that by 
taking E? of this form, the four diffeomorphism constraints (9.27)-(9.28) are automatically 
satisfied. 
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Substituting (9.34) into the vector constraint (9.28), we get 

E\F ah = i{r% 3 B h ^ hc B™ = 0, (9.35) 

where we have used the fact that (ip -1 )^ is symmetric in % and j while is anti-symmetric 
in b and c. Similarly, substituting (9.34) into the scalar constraint (9.27), we get 

e ijk^b Fabk = e ^E?E b jrbbc B c k 

= -te^E^E^^E? 

= -iqe^e^l (9.36) 

= -2un# 
= o, 

where we have used the fact that ipij is trace-free. Thus, the four diffeomorphism constraints 
(9.27)-(9.28) are automatically satisfied for Ef having the form given by (9.34). That this 
is the most general solution follows if Bf is non-degenerate. Then for a given A l a , Ef will 
have 5 degrees of freedom (per space point) corresponding to the 5 degrees of freedom of the 
symmetric trace-free tensor ipij. 

What remains to be solved is the Gauss constraint (9.29), which in terms of B m and 
can be written as 

= V a Et 

= tV a ((4>-%B a i) (9.37) 
= iB a iV a (i>-%. 

To obtain the last line of (9.37), we used the Bianchi identity V a B aj = \rj abc V [a F 3 bc] = 0. 
10. Discussion 

Let me begin by briefly summarizing the main results reviewed in this paper. 

I. The standard Einstein-Hilbert theory is a geometrodynamical theory of gravity in 
which a spacetime metric is the fundamental field variable. The phase space variables 
consist of a positive-definite metric q a b and its canonically conjugate momentum p ab . 
These variables are subject to a set of 1st class constraints, which are non-polynomial 
in q a t, and which have a Poisson bracket algebra involving structure functions. This 
theory is valid in n + 1 dimensions. 
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2. The 2+1 Palatini theory is a connection dynamical theory defined for any Lie group 
G. The fundamental field variables consist of a £<3-valued connection 1-form and a 
£g-valued covector field. The phase space is coordinatized by a connection 1-form 

and its canonically conjugate momentum (or "electric field") Ef. These are fields 
defined on a 2-manifold £, and they are subject to a set of 1st class constraints. The 
constraints are polynomial in (A T a , Ef) and provide a representation of the Lie algebra 
of the inhomogeneous Lie group associated with G. One recovers 2+1 gravity by taking 
G = 50(2,1). 

3. Chern-Simons theory is a connection dynamical theory defined for any Lie group that 
admits an invariant, non-degenerate bilinear form. In 2+1 dimensions, the fundamental 
field variable is a Lie algebra-valued connection 1-form, and the phase space variables 
are A\ — the pull-back of the field variable to the 2-dimensional hypersurface S. There 
are 1st class constraints, which are polynomial in A\ and provide a representation of 
the defining Lie algebra. Chern-Simons theory is related to 2+1 Palatini theories as 
follows: (i) 2+1 Palatini theory based on any Lie group G is equivalent to Chern- 
Simons theory based on the inhomogeneous Lie group IG associated with G; and (ii) 
the reduced phase space of Chern-Simons theory based on a Lie group G is the same 
as the reduced configuration space of the 2+1 Palatini theory based on the same G. 
As a special case of (i), 2+1 gravity is equivalent to Chern-Simons theory based on 
the 2+1 dimensional Poincare group I SO (2, 1). 

4. One can couple matter to 2+1 gravity via the 2+1 Palatini action. 2+1 Palatini theory 
coupled to a cosmological constant A is defined for any Lie group G that admits 
an invariant, totally antisymmetric tensor e IJK . This theory is equivalent to 2+1 
dimensional Chern-Simons theory based on the A-deformation of G. As a special case, 
2+1 gravity coupled to a cosmological constant is equivalent to Chern-Simons theory 
based on 50(3, 1) or 50(2, 2) (depending on the sign of A). 2+1 Palatini theory can 
also be coupled to matter fields with local degrees of freedom provided G = 50(2, 1). 
The constraints remain polynomial in the canonically conjugate variables and form a 
1st class set. However, due to the presence of structure functions, they no longer form 
a Lie algebra. 

5. The 3+1 Palatini theory is a geometrodynamical theory of 3+1 gravity in which a 
co-tetrad and a Lorentz connection 1-form are the fundamental field variables. Due 
to the presence of 2nd class constraints, the Hamiltonian formulation of this theory 
reduces to that of the standard Einstein-Hilbert theory in 3+1 dimensions. Unlike the 
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2+1 Palatini theory, the 3+1 Palatini theory does not provide a connection dynamical 
theory of 3+1 gravity. 

6. The self-dual theory is a connection dynamical theory of complex 3+1 gravity in which 
a complex co-tetrad and a self-dual connection 1-form are the fundamental field vari- 
ables. The phase space variables consist of an CC S o(3)~ valued connection 1-form A\ 
and its canonically conjugate momentum (or "electric field") both defined on a 
3-manifold E. These variables are subject to a set of 1st class constraints, which are 
polynomial in (A l a , Ef) but which have a Poisson bracket algebra involving structure 
functions. In a solution, Ef is a (densitized) spatial triad. Since none of the equations 
involve the inverse of Ef, the self-dual theory makes sense even if Ef is non-invertible. 
Thus, the self-dual theory provides an extension of complex general relativity that is 
valid even when the induced spatial metric q ab (= EfE b% ) is degenerate. One must 
impose reality conditions to recover real general relativity. 

7. One can couple matter to complex 3+1 gravity via the self-dual action. The constraints 
remain polynomial in the canonically conjugate variables and form a 1st class set. Since 
none of the equations involves the inverse of Ef, the self-dual theory coupled to matter 
provides an extension of complex general relativity coupled to matter that includes 
degenerate spatial metrics. Reality conditions must be imposed to recover the real 
theory. 

8. The pure spin-connection formulation of general relativity is a connection dynamical 
theory of complex 3+1 gravity in which a CCso (3) -valued connection 1-form and a 
scalar density of weight —1 are the fundamental field variables. The Hamiltonian 
formulation of this theory is equivalent to that of the self-dual theory provided the self- 
dual part of the Weyl tensor is non-degenerate. When the "magnetic" field associated 
with the connection 1-form A l a is non-degenerate, one can write down the most general 
solution to the four diffeomorphism constraints. Reality conditions must be imposed 
to recover the real theory. 

So what can we conclude from all these results? Is gravity a theory of a metric or a 
connection? In other words, is gravity a theory of geometry, where the fundamental variable 
is a spacetime metric which specifies distances between nearby events, or is it a theory of 
curvature, where the fundamental variable is a connection 1-form which tells us how to 
parallel propagate vectors around closed loops? The answer: Either. As far as the classical 
equations of motion are concerned, both a metric and a connection describe gravity equally 
well in 2+1 and 3+1 dimensions. Neither metric nor connection is preferred. 
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As we have shown in this review and have summarized above, 2+1 and 3+1 gravity admit 
formulations in terms of metrics and connections. But despite the apparent differences (i.e., 
the different actions and field variables; the different Hamiltonian formulations and canon- 
ically conjugate momenta; and the possiblity of extending the theories to include arbitrary 
gauge groups and solutions with degenerate spatial metrics), we have seen that the classical 
equations of motion for all these formulations are the same. For instance, we saw that the 
2+1 Palatini theory reproduces vacuum 2+1 gravity when we choose G = 50(2, 1) and solve 
the equation of motion for the connection. Similarly, we saw that Chern-Simons theory re- 
produces 2+1 gravity coupled to a cosmological constant A (> 0) when we choose the gauge 
group to be 50(2,2). At the level of field equations, all the theories are mathematically 
equivalent. The difference between the theories is, instead, one of emphasis. 

Now such a small change may not seem, at first, to be worth all the effort. Recall that 
the shift in emphasis from metric to connection came only after we successively analyzed the 
Einstein-Hilbert, Palatini, and Chern-Simons theories in 2+1 dimensions, and the Einstein- 
Hilbert, Palatini, self-dual, and pure spin-connection theories in 3+1 dimesions. This analysis 
required a fair amount of work and, as we argued in the previous paragraph, did not lead 
to anything particularly new at the classical level modulo, of course, the extensions of the 
theories to include arbitrary gauge groups and solutions with degenerate spatial metrics. 
But as soon as we turn to quantum theory and consider the recent results that have been 
obtained there, the question as to whether the shift in emphasis from metric to connection 
was worth the effort has a simple affirmative answer. Yes! Indeed, almost all of the recent 
advances in quantum general relativity can be traced back to this change of emphasis. As 
mentioned in the introduction, Witten H used the equivalence of the 2+1 Palatini theory 
based on 50(2, 1) with Chern-Simons theory based on 750(2, 1) to quantize 2+1 gravity. 
Others (e.g., Carlip Jj4], |25j and Anderson ||26|| ) are now using Witten's quantization to 
analyze the problem of time in the 2+1 theory. In 3+1 dimensions, Jacobson, Rovelli, and 
Smolin || |7|] took advantage of the simplicity of the constraint equations in the self-dual 
formulation of 3+1 gravity to solve the quantum constraints exactly — something that nobody 
could accomplish for the quantum version of the scalar constraint in the traditional metric 
variables. And the list goes on. (See, e.g., @, § and JL7|, [Tj| [L^, |20| for more details.) Where 
this list will end, and whether or not the change in emphasis from metric to connection 
will lead to a mathematically consistent and physically reasonable quantum theory of 3+1 
gravity, remains to be seen. 
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